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THE THEORY OF MICRO-METEORITES.* PART I. IN AN 
ISOTHERMAL ATMOSPHERE 
By Frep L. 
HARVARD COLLEGE OBSERVATORY 
Communicated by Harlow Shapley, October 6, 1950 


Introduction.In 1946, during the great Giacobinid meteor shower, 
H. E. Landsberg! collected several small magnetic particles that apparently 
were associated with the shower. Since some of these particles, a few 
microns in length, were extremely angular in shape (wedge-shaped and 
opaque), it seemed unlikely that they could have been the end-products 
of vaporizing meteors. Landsberg concluded that they must have been 
stopped by the atmosphere without being heated above their melting-points. 
As a result of his suggestion I have developed the present theory to in- 
vestigate the process whereby temperature radiation can dissipate the 
energy gained by encounters with atmospheric molecules sufficiently 
rapidly to permit finite meteoric particles to be stopped without melting. 
Some basic concepts of this theory have been discussed by E. Opik? and an 
application made in the case of an isothermal atmosphere. 

The term micro-meteorite appears to be an appropriate designation for one 
of these small particles. 

In every sense the micro-meteorites represent the lower extreme to the 
ascending sequence embracing meteor, fireball and meteoritic crater formation. 
Hence, the theory is a limited meteor theory, partially applicable to the 
unobservable beginning of a meteor. 

We may now assume and later (Part IT) prove that interaction between 
the air molecules striking and those leaving the micro-meteorite may be 
neglected. The molecular mean free paths, even after correction for the 
relatively slow velocity of air molecules thermally emitted, are greater than 
the linear dimensions of the micro-meteorite. 

Let us suppose that the micro-meteorite presents a certain surface 
area, A, to the atmosphere, which it encounters with velocity, V. This 
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surface may be an average frontal area in case the body is rotating, or a 
fixed area in case the body does not rotate. If A(@, 6), a function of the 
ordinary spherical coordinates, represents the actual cross-sectional profile 
of the micro-meteorite as seen from direction ¢ and @, then the average 
frontal area, A, is given by 


j 
A, = / / Alo, 4) sin de dé. (1) 


This average area, A,, may or may not equal the frontal area A. The 
thermal radiating area, B, however, will be specific for a particle of a given 
shape; it is given by 


B 4A). (2) 


Let us assume that the temperature of the micro-meteorite is at all times 
uniform over the area B. This assumption is equivalent to an assumption 
that the heat conductivity is infinite or that the heat capacity is zero. 
We shall first develop the theory on the assumption that the heat capacity 
is negligible and later investigate the nature of the error made. 

The air molecules impinge on the forward surface with relative velocity, 
l’, because, by definition, the mean free path of the outgoing molecules 
relatively to the moving body is larger than its linear dimensions. To 
evaluate the energy transfer to the surface of the body we may make use 
of the concept of the accommodation coefficient, a, which is defined in 
terms of the kinetic energy of the air molecules, referred to the coordinate 
system of the moving body. The accommodation coefficient is, then, the 
actual loss of kinetic energy by the air molecules, as a result of the encoun- 
ter, divided by the loss if all of them were momentarily to adhere to the sur 
face and be re-emitted at the thermal velocity corresponding to the surface - 
temperature. Since the thermal energies, both original and at the surface 
temperature, are relatively small compared to the kinetic energy at ve- 
locity, V’, these energies may be neglected with an error less than 1% 
(minimum Vo = 11.2 km.‘see.). Hence a@ represents the fraction of the 
molecular energy at velocity, V, that is transmitted to the micro-meteorite. 

Since the air molecules will encounter the body with relative energies 
of the order of 8 to SOO electron volts, while the work function of the sur- 
face will be only a few volts, the molecules will certainly penetrate the sur- 
face for several molecular layers except at the lowest velocities. We must 
conclude that few of them will leave with high velocities; the losses by dis 
sociation, excitation and ionization can be only a few volts. Hence a 
must be nearly unity at most velocities 

In air of density, p, the micro-meteorite will meet in time, df, an air 


mass, given by 
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dm, = ApV dt. 


The corresponding energy gain, dE,, amounts to 


dE, = dm, = ~ ApV? dt. (4) 


Part of this energy will be utilized in raising the temperature of the 
meteoroid, part radiated by black-(or gray)-body radiation, part used in 
dissociation, excitation and ionization and perhaps part used in disengaging 
material from the surface. If the accommodation coefficient is defined to 
include the dissociation, excitation and ionization and if vaporization is 
negligible, we may deal here explicitly with only the heating and radiation 
terms. 

We may assume that the meteoroid was previously in temperature 
equilibrium with the night side of the Earth (or Sun and Earth during 
the day) at temperature 7). With a gray-body emissivity coefficient of 
8, the loss of energy by radiation, dE,, of a surface at temperature T, is 


dE, = BBa (T,* — To')dt, (5) 


where @ is the Stefan-Boltzmann constant. 

If the meteoroid is small (s < 1 cm.), of mass, m, and if the coefficient 
of heat conductivity is at all comparable to that of ordinary rocks, the 
internal temperature should differ negligibly from the surface temperature 
in time intervals somewhat smaller than one second. The permissible limits 
to this assumption will be discussed later. If, then, the heat capacity per 
gram is C,, the temperature will vary as , 


mCdT, = dE, dE,. (6) 
By equations (4) and (5), equation (6) becomes 


mC, = ApV® — BBo(T,‘ — T>'). (7) 

The precise conditions under which the heat capacity in the left member 
of equation (7) can be neglected are not apparent a priori. We can, how- 
ever, easily determine a resultant rough limit to the dimensions of the 
micro-meteorite and later study the question more thoroughly. The maxi- 
mum temperature, 7,,, to which a micro-meteorite can be heated without 
appreciable vaporization is just below the melting point of the least re- 
fractory material in the meteoroid, approximately 1200°K. to 1700°K. for 
typical stones.* Iron, iron oxides and silica also fall within this range. 
The temperature rise from the equilibrium temperature at the Earth to 
T,, is relatively large. GeneraHy this rise will occur over a considerable 
distance through the atmosphere since the radiation varies as 7,‘ while the 
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heating is proportional to atmospheric density. We will find that the 
heat capacity is negligible for fast micro-meteorites because they are so 
small. For slow micro-meteorites an atmospheric density increase of a 
factor of two will require roughly a half-second of atmospheric traverse at 
normal incidence. 

Hence, if the heat radiation in an interval of the order of half-second near 
maximum temperature is large compared to the heat required to raise the- 
temperature from 7) to T,, we may safely neglect the heat capacity in 
equation (7). The condition just defined is 


The limiting radius, s, for a spherical meteor of density, p,, is, from 
equation (8) 


+ + Tale? + T,*). (9) 


The right member of equation (9) is of the order of 0.01 cm. for an iron 
or stony meteorite. Hence we may safely ignore the heat capacity of 
micro-meteorites of radii less than 10 microns. 

With our current assumptions then, we may set the left-hand member 
of equation (7) equal to zero and determine the surface temperature of a 
micro-meteorite as a function of its velocity and the atmospheric density. 
The result is 

aApV? 


T 
28Be 


(10) 

While the temperature of the micro-meteorite is rising, its velocity is 
being reduced by atmospheric resistance. We may define the drag co- 
efficient, D, by the equation 


ADpV? 


2 


md V = (11) 

The drag coefficient will include a major component, of the order of 
2 numerically, if we assume that air particles momentarily adhere to the 
meteoroid and are reemitted isotropically. A smaller additional term 
will arise from the non-isotropic reemission according to the hypothesis of 
Tsien‘ and Miss Heineman.* No term will arise from evaporation of ma- 
terial, however, since we assume (perhaps erroneously at high velocities) 
that evaporation is negligible. 

A rigorous inclusion of the reemission drag term would involve the sur- 
face temperature of the micro-meteorite. Since the effect is rather small, 
we may include it approximately after a*solution has been made with D 
constant, 
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The simultaneous solution of equations (10) and (11) will provide a re- 
lation among the temperature, velocity and time, if the atmospheric 
density can be expressed as a function of the time. A complicated, vari- 
able and poorly known relation between atmospheric density and height 
actually exists at the atmospheric heights concerned. In the following 
sections the fundamental equations will be solved for certain types of such 
relations. 

The Solution in an Isothermal Atmosphere.—Let us assume that the at- 
mosphere is isothermal and that the mean molecular weight of the air is 
also constant. Then the atmospheric density, o, varies with the height, 
h, above sea level according to the relation 


p= pe, (12) 
where 4, the logarithmic density gradient, is a positive constant and pp 
is constant. 

Let us neglect the curvature of the Earth. Then for a particle entering 
the atmosphere with velocity, V, from an apparent radiant at zenith dis- 
tance, Z, the time and height are related by the equation 

dh = —V cos Z dt. (13) 


Since meteoric velocities are so large we may neglect the effect of the 
Earth's gravity, assuming that the pertinent part of the trajectory is a 
straight line, and that the velocity is unaffected by gravity.* 

Equations (11), (12) and (13) then lead to the following differential 
relation between velocity and height: 


(14) 


We have already assumed that all of the quantities in equation (14) ex- 
cept V and h may be treated as constants. Hence we may integrate equa- 
tion (14) between the limits of V,, the initial velocity at great heights 
essentially infinity, to any pertinent velocity, V, and height, 4. The in- 
tegral of equation (14), when combined with equation (12) to eliminate 
the height explicitly, yields the following relation between the velocity 
and atmospheric density: 


ADp 


log (V/V...) = — (15) 


2bm cos Z 


Now the micro-meteorite can traverse the atmosphere undamaged only 
if the velocity is reduced sufficiently to prevent the surface temperature 
from exceeding a critical temperature, 7, somewhat below the melting 
point of the meteoritic materials, This restriction can be imposed if we 
eliminate the density, p, between Equations (15) and (10) and then impose 
a maxima) condition on T,. The first step gives 
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we alm cos Z : 
~ bog (V/V,). (16) 
To determine the maximum temperature from equation (16) we must 
make assumptions as to the functional forms of the accommodation co- 
efficient, a, and the emissivity factor, 8. In view of our lack of knowl- 
edge of the detailed physical and chemical structure of the micro-meteorites 
and, even with this knowledge, the uncertainty in the physical laws for a 
and 8, we may as well adopt constant values for these quantities as well as 
for D. Hence, from the derivative of equation (16) the maximum tem 
perature is reached at the critical velocity, V,, given by 


log V.) = (17a) 


V. = = 0.7165 (17b) 


where ¢ is the base of the natural logarithms. 

We may now solve for the ratio of the mass to the effective surface 
area, m/B, of the micro-meteorite from equation (16) for the maximum 
temperature, 7. The resulting critical ratio 1s 

seBoD(T,,* — 
m/B = (18a) 
ab cos 

Equation (18a) represents the maximum value of m, B for a micro-mete- 
orite that is to be stopped, undamaged by the atmosphere. The maxi- 
mum radius for a spherical particle of radius, s (sphere) and density, p,, 1s 
then 


— 
{ LSb) 
cos ZV,,* 


s (sphere) = 


a result corresponding closely to Opik's Equation (53).* 

It is of interest to note from equations (18) that the maximum particle 
dimension depends directly upon the critical temperature to the fourth 
power and upon the inverse cube of the original velocity. Low velocity and 
high melting point strongly favor the passage of such a particle through 
the atmosphere. Irregular or elongated shape favors the process for a 
given meteoric mass. Because of the cos Z term, larger masses may enter 
at lower angles of incidence. 

For a rapidly spinning micro-meteorite in the shape of a right circular 
cylinder of length / and radius s (cyl.), the maximum radius in terms of the 
comparable sphere (equation 1Sb) of the same density and m/B ratio is 
given by 
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2/s + 
s (cyl.) = s (sphere) X ;( ; ‘) (cyl.). (18e) 


If the rapidly rotating cylindrical meteorite is to have the same radius as 
the spherical one, then the length must be equal to the diameter of the 
sphere. An extremely long cylindrical meteorite (/ > s) will be stopped 
when its radius is 2/3 or less the corresponding spherical radius (equation 
18b). The maximum permissible dimensions of micro-meteorites with 
other shapes or orientations can be determined by means of equations (1), 
(2) and (18). 

In applying equations (18) approximately to an atmosphere of variable 
temperature it is necessary to obtain an appropriate value of the logarithmic 
density gradient, 6. The critical altitude is at the point of maximum tem- 
perature for the micro-meteorite. From equations (15) and (17) we find 
the critical density, p,,, at maximum temperature T,, as follows: 

2bm cos Z 


If the frontal area of the meteorite, A, is approximated by the average 
area, A;, equations (2) and (18a) transform equation (19) into 


SeB8o(T,,* 
aV.? 

The value of 6 and the height corresponding to p,, may be derived from 
some standard atmosphere. The fact that the shape and mass of the 
micro-meteorite and cos Z do not enter equation (20) is rather surprising. 
In case there is reason to believe that the actual frontal area differs from 
the average cross-sectional area, a correction term can easily be applied 
for this factor in equation (20). 

The drag coefficient, D, is a vital factor in equations (18) and (19). - In 
case the air molecules impinging on the surface of the micro-meteorite 
are momentarily “captured” by the surface, but very quickly reemitted 
with energies corresponding to the temperature of the surface, D is equal 
to 2 with a small correction term for the Tsien-Heineman effect mentioned 
above. This situation corresponds to that when the accommodation co- 
efficient, a, is exactly unity. M. L. Wiedman’ finds that the accommoda- 
tion coefficient for air on various metals is generally close to 0.9. At the 
higher energies here considered, there is little doubt that the air molecules 
will largely penetrate the surface molecular layers and be momentarily 
captured. Since only a few can combine chemically and since condensa- 
tion is impossible at the temperatures considered, there appears to be 
little question that essentially all of the air molecules will be captured and 
then reemitted with an average velocity, v7, corresponding to the tem- 
perature of the surface. The average velocity is given by 


Pm ( 19) 


Pm (20) 
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(vy)? = (21) 
where &, is the gas constant and yw is the molecular weight of the gas. 

If the micro-meteorite ss assumed to be spherical and the reemission 
is uniform with respect to solid angle from a small area of the surface, 
the momentum transfer to a hemisphere per average molecule is 4yrr/9. 
Since the time lag of reemission must be small compared to a possible rota- 
tion period of the micro-meteorite, this momentum transfer will oppose 
the motion of the body. If we further accept D = 2 as the best approxi- 
mation to the drag coefficient without including reemission, the value 
of D becomes: 


D = 2[1 + (22 
A numerical calculation shows that the tp-correction in equation (22) 
amounts at most to a few per cent in the velocity and temperature ranges 
under consideration, Most of the deceleration of interest occurs fairly 
near the critical region of maximum temperature. In adopting a con- 
stant value of D, therefore, we may use vy applying at the maximum tem- 
perature and compensate this error somewhat by adopting V = V,, in 


equation (22). 

It would be desirable to make some correction in the drag coefficient 
for the fact that some of the impinging air molecules must deviate from 
the assumed process of penetration and thermal reemission, but we have no 
detailed information as to the surface characteristics of micro-meteorites. 
The relative error made in D by this omission is less than 1 — a. 

The remaining undetermined constant in Equation (18) is the emissivity 
factor, 8. Nothing precise can be known about this factor until micro- 
meteorites have been studied carefully in the laboratory. Probably § is 
very near unity. For want of better information we may set a = 6 = 1, 
so that their effects cancel in equation (18). We may adopt 7,, = 1600°K., 
and 7, ~ 300°K. The remaining quantities to be specified in equation 
(18a) concern the atmosphere. The Tentative Standard Atmosphere 
of the National Advisory Committee for Aeronautics*® represents good 
modern estimates of upper-atmospheric densities and will be used for the 
present calculation 

For a velocity of 23 km. sec., corresponding to the 1947 Giacobinid 
Meteor Shower, the atmospheric density at which a micro-meteorite of 
limiting dimensions attains maximum temperature = 8.2 
gm. em.") occurs near the 112-km. level, by application of equation (20). 
With an air temperature of 346°K. and b = 1.1 K 10>-* cm.~! the maxi- 
mum radius of spherical iron particles for cos Z = 0.45 are then calculated 
to be about 4 microns. Long cylindrical particles might have a diameter 
as great as 6 microns, in good agreement with Landsberg’s observations 


mentioned earlier 
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Part II of this paper will concern the solution of limiting dimensions 
for micro-meteorites in an atmosphere of constant temperature gradient 
and in the general case. A critical discussion of the assumptions adopted 
here will be included as well as some elaboration of the numerical results. 

* This paper was written as a part of the investigation for United States Naval Ord- 
nance Contract-10449-05512. 

' Pop. Ast., $5, 322 (1947). 

* Pub. Unw. Tartu, 29, No. 5, 51 (1937). 

* See Daly, R. A., Igneous Rocks and the Depths of the Earth, McGraw-Hill Book Co., 
1933, p. 65. 

* Tsien, Hsue-Shen, J. Aero. Sct., 13, 653 (1946). 

* Heineman, M., Comm. A pl. Math., 1, 259 (1948). 

* See, e. g.; Whipple, F. L., Proc. Am. Phil. Soc., 79, 499 (1938). 

? See Tsien, loc. cit. 

* Warfield, C. N.,. NACA, Tech. Note, No. 1200, Langley Field, Jan., 1947. 


MERCAPTAN-INDUCED COAGULATION OF PROTEINS* 
By CHARLES HucGins, DONALD FRASER TAPLEY AND ELwoop V. JENSEN 


From THE NATHAN GOLDBLATT Memorial Hosprral FOR NBOPLASTIC DIsRASES OF THE 
Untrversity or Curcaco, Curcaco, IL. 


Communicated, October 29, 1950 


In this paper it will be demonstrated that mercaptans possess the 
property of inducing coagulation of certain proteins at room temperature 
and neutrality. This effect throws light on the importance of —-S—S— 
bonds in the intramolecular folding of several proteins in the native state. 

While the dispersive action of mercaptans on keratin is well known, 
the coagulative action on soluble proteins apparently has not been de- 
scribed. Goddard and Michaelis’ found that wool readily dissolves at 
room temperature in thioglycollate solutions at pH 10-13 because of 
reduction of disulfide bonds. The following reaction'~* occurs in the thiol- 
disulfide system: 

RS—SR + 2R’SH 2RSH + R’'S—SR’ 
Since greater alkalinity was required in the experiments of Goddard and 
Michaelis than was needed for the reduction of cystine by thiols, they" 
postulated that an additional intramolecular bridge must be broken before 
disulfides of keratin can be reduced; they assumed that the second bridge 
was salt-like in character and that it was disrupted by removal ef a proton 
from the amino group in alkaline solution. Jones and Mecham’ observed 
that a variety of keratins were dispersed at pH 7 and 40°C., by adding 
urea or anionic detergents to the mercaptan solution. Mirsky and Anson? 
found that thioglycollate in excess completely reduces disulfide groups in 
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denatured plasma albumin and egg albumin; in contrast sodium sulfite 
reduces only one-half of the cystine present. 

Method.—Crystalline bovine plasma albumin and egg albumin and 
refined bovine y-globulin (95 per cent homogeneous on electrophoresis) 
were studied. Acetylated bovine albumin was prepared by the method of 
Hughes as described by Olcott and Fraenkel-Conrat.* Analysis indicated 
that 78 per cent of the free amino groups were acetylated. 

Unless otherwise stated the proteins were dissolved in 0.1 M phosphate 
buffer, pH 7.4, and the experiments were carried out in 13 XK 100 mm. 


ATE 


PIGURE 1 FIGURE 2 


Rates of development of turbidity Differential rates of development of 


in plasma albumin induced by 0.2 M 
solutions of monothioglycol, and sodium 
thioglycollate. In each case the system 
consisted of bovine plasma albumin 50 
mg., Mercaptan compound 2 m. eq. and 
0.1 M phosphate buffer, pH 7.4, to make 
the final volume 10 ml.; temperature 
37° 


glass tubes at 37°. 


eter at A = 660 miullimicrons, using 19- X 150-mm. tubes. 


turbidity in various refined protein 
fractions induced by 0.2 M mercapto- 
ethanol. The system consisted in 
each case of 50 mg. protein, mercapto- 
ethanol 2 m. eq., and 0.1 M phosphate 
buffer pH 7.4 to make the final volume 
10 ml.; temperature 37°. 


Turbidity was determined in a Coleman spectrophotom- 


A protein 


solution was considered to be coagulated’ when it failed to pass through 

a clean dry stainless steel screen (8 mesh, wire diameter 0.7 mm.). 
Experimenial.-When a mercaptan is added to a solution of bovine 

plasma albumin turbidity soon develops and increases to a maximum 


(Fig. 1), then decreasing shghtly. 
until a firm opaque gel forms. 


The viscosity of the solution increases 
Egg albumin and y-globulin behave in a 


sitnilar fashion except that the reaction is much slower (Fig. 2). Dialysis 
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for 5 days against water produces no apparent change in these white 
mercaptan-induced gels. 

An inverse relationship exists between the protein concentration and the 
minimal quantity of thioglycollate required to form a gel; at pH 7.4 and 
37° the thioglycollate requirement decreased progressively from a con- 
centration of 0.34 M to 0.08 M as the albumin concentration was increased 
from 1.2to5percent. At 5°, however, neither turbidity nor gels developed 
in 5 per cent albumin even after 72 hr. in 0.4 M thioglycollate. 

In solutions of 7 per cent albumin, turbidity with subsequent gelation 
developed with the following agents: mercaptoethanol (monothioglycol), 
sodium thioglycollate, reduced glutathione, cysteine and sodiun hydro- 
sulfite. Mercaptoethanol produced turbidity more rapidly (Fig. 1) and 
gel formation at a lower concentration than did the other reducing agents 
employed. No gels or turbidity occurred in 7 per cent albumin solutions 
with the following reagents in concentration up to 0.4 M: sodium cyanide, 
potassium ferrocyanide, hydroquinone or sodium sulfite. Slight turbidity 
but no gel formation occurred in 0.4 M sodium ascorbate. 

Coagulation of albumin by thiol compounds is favored by an increase 
in pH of the mixture. With increasing pH the minimum concentration of 
sodium thioglycollate required to form a gel im 2.5 per cent solutions of 
plasma albumin decreased progressively from 0.18 M at pH 7.0 to 0.06 M 
at pH 8.4. No monophasic gels were observed in albumin solutions treated 
with thioglycollic acid at pH values below 6.8, Between pH 5 and 6.6 
the protein solution remained clear in 0.2 M thioglycollic acid while at 
pH 4.5 a synerizing precipitate formed. 

Variation of the salt concentration over the range 0.1 to 1.1 Af had no 
apparent effect on mercaptan-induced gelation of albumin. At salt con- 
centrations less than 0.02 M gel formation was markedly inhibited, while 
at salt concentrations of 2-3 M the coagula formed were quite friable. 
Oxygen apparently plays no part in the mercaptan-induced coagulation of 
albumin; solutions degassed by repeated freezing and thawing in vacuo 
and sealed in the limbs of A-tubes evacuated to 107-* mm. of Hg before 
mixing behaved identically with control mixtures exposed to air. 

Partial acetylation (78 per cent) of the amino groups of plasma albumin 
resulted in a twofold increase in the mercaptan concentration required to 
form a gel; in a 2.5 per cent protein solution at pH 7.4, 0.27 M thioglycol- 
late was required instead of 0.15 M as with unmodified albumin. The gels 
thus formed from acetylated albumin were markedly less turbid than those 
formed from unmodified albumin itself. 

The opaque gels formed by the mercaptans are rapidly dissolved by the 
addition of urea in final concentration of 5 M or greater with nearly com- 
plete disappearance of the turbidity. Now urea itself possesses the 
ability to induce gelation of albumins and certain other proteins.» * These 
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coagula are quite different from mercaptan imduced gels, being translucent 
and very cohesive. The protein gels induced by urea are dissolved rapidly 
by the addition of mercaptans.” 

Discussion. Two types of chemical bonds have been established as 
forces which maintain the protein molecules in specific configuration, 
namely, hydrogen bonds and disulfide bridges. In the theory of Mirsky 
and Pauling'' the polypeptide chain of a native protein is folded and re- 
tained in its unique configuration by hydrogen bonding between the 
peptide nitrogen and oxygen atoms and also between the free amino and 
carboxyl groups of the diamino and dicarboxylic amino acid residues. 
Other workers have postulated that cystine’* must form one of the side 
links between adjacent polypeptide chains in the protein molecule. 

Classes of compounds which disrupt either type of linkage lead to 
the denaturation of a protein such as plasma albumin. Amides are 
well-known hydrogen-bond-forming substances which can rupture the 
hydrogen bonds involved in maintaining the protein in the native state; 
mercaptans reduce —S-—S-—- linkages. With each type of reagent the 
denatured molecules aggregate, the gels varying in properties because of 
differences in the denatured protein molecule. 

The coagulation of albumin by urea has recently been interpreted" as 
follows: Concentrated amides by competition with intramolecular hydro- 
gen bonds of the native protein disrupt its structure and make disulfide 
groups reactive to small concentrations of sulfhydryl. By means of a 
chain reaction between protein —-SH and —-S-—-S— groups a reticulum is 
then knitted together consisting of intermolecular protein disulfide bonds 
which hold together the unfolded protein chains. Solvent is bound to the 
protein within this framework so that clear firm gels result. 

Mercaptans appear to coagulate certain soluble proteins by a different 
mechanism, namely, disruption of the protein configuration through the 
cleavage of disulfide bonds followed by aggregation of the denatured protein 
through hydrogen bonds. Plasma albumin undergoes these changes more 
rapidly than y-globulin or egg albumin implying that the disulfide groups 
of albumin are more readily available to reduction of —-SH groups than in 
the other proteins. 

The ready dispersion of the opaque mercaptan-induced gels upon the 
addition of urea indicates that aggregation in these gels is chiefly through 
intermolecular hydrogen bonding. Further evidence is derived from the 
deficient capacity of acetylated albumin to coagulate in the presence of 
mercaptans, due to a decrease in the number of potential hydrogen bond 
formers, the amino groups. The absence of any oxygen effect in mer- 
captan-induced gelation makes it doubtful that reoxidation of reduced 
disulfide bonds plays a major part in the aggregation. 

Summary.—Mercaptans and hydrosulfite denature and subsequently 
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coagulate bovine plasma albumin, y-globulin and egg albumin to form 
opaque gels at room temperature and neutrality. The reaction with 
bovine plasma albumin is many times more rapid than with the other two 
proteins. Gels so formed rapidly dissolve in concentrated solutions of 
amides such as urea. The capacity of albumin to form this type of gel is 
decreased by partial acetylation of the albumin. The evidence indicates 
that mercaptan-induced gelation occurs through hydrogen bonding be- 
tween the extended denatured protein molecules. 


* This investigation was aided by grants from Mr. Ben May, Mobile, Ala., from the 
Daisy Schwimmer Fund, from Mr. Ernest G. Shinner, from the Anna S, Goldberg 
Scholarship Fund and from the American Cancer Society recommended by the Com- 
mittee on Growth of the National Research Council. 
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1HE ACTION SPECTRUM OF THE INHIBITION OF GALACTO- 
ZYMASE PRODUCTION BY ULTRAVIOLET LIGHT* 


By Paut A. Swensont 
DEPARTMENT OF STANFORD University, STANFORD, CALIFORNIA 
Communicated by V. C, Twitty, October 16, 1950 


Certain yeasts, when placed in contact with galactose, acquire the 
ability to use that substrate after a period of time. Such a response is 
called enzymatic adaptation and will take place even in the absence of a 
nitrogen source. In recent years, with the realization that at least one of 
the components of any enzyme is a protein, it has been recognized that in 
adaptive enzyme formation we have a case where specific protein material 
is synthesized by resting cells. 

It has been demonstrated that ultra-violet radiations from a sterilamp 
inhibit the aerobic adaptation of yeast to galactose.'’ Determination of 
an action spectrum for this effect might enable us to identify the cellular 
constituent affected by the radiation. The problem is of theoretical in- 


I. 


4 
A 
7 
4 
| 
ae 
} 
i 
ae 
‘ 


700 BIOCHEMISTRY; P. A. SWENSON Proc. N.A.S 


terest because a clue might be found to the mechanism of adaptive enzyme 
synthesis. An action spectrum has been determined and is reported herein. 

The theoretical basis and the limitations of the action spectrum method 
are discussed by Loofbourow.’ The reciprocals of the energy required to 
produce a given effect when plotted against the wave-length give a curve 
which is called the action spectrum, This is thought to correspond to the 
absorption spectrum of some photolabile compound responsible for the 
biological effect being considered; the absorption spectra of known com- 
pounds are matched with the action spectrum. A summary of ultra- 
violet action spectra is given by Giese.* 

Experimental.—In the present experiments the same strain of yeast was 
used as in the previous experiments.' Cultures were grown at 28°C. for 
24 hours in liquid medium (30 ml. 3°) dextrose, 1% yeast extract in 250 
ml. flasks) on a mechanical shaker. The cells, washed as before, were sus- 
pended for experiments in M/20 KH, PO, in a concentration of 2 X 10° 
cells per ml. 

TABLE 1 
MINIMAL Doses Regurrep 10 Prevenr GALACTOZYMASE FORMATION FoR AT LEAST 
Four Hours 


Wave 
ceworm, A 


2345 3374 
2537 1687 
2654 1687 
280M 3374 
3025 13,596 
3130 > 54,384 
Ultra-violet light obtained from a quartz mercury are was passed through 
a quartz monochromator, and the intensity of the light was measured with 
a thermopile galvanometer system as described elsewhere.‘ The mono- 
chromatic light struck the suspension of cells from below. The volume 
of cells irradiated was 0.06 ml. Stirring during irradiation was performed 
by a moving jet of compressed air saturated with water vapor. Pre- 
cautions were taken against evaporation; if evaporation was excessive, 
the sample was discarded. Once irradiation had begun, the sample was 
handled only in red or yellow light to prevent photoreactivation. 
Following irradiation, the sample was transferred to the reaction cup of 
a capillary microrespirometer of a modified Kirk-Stern type,’ developed 
by the author and Mr. Robert Perthel of the Stanford Research Institute. 
A 0.06-ml. sample of 4°) galactose solution was added to the suspension 
with a braking pipette, and the course of the respiration followed by re- 
cording the progress of a kerosene drop in the capillary tube. Thermal 
regulation of the system was accomplished by submerging the cups at 
either end of the U-shaped capillary in a vessel of water at room tempera- 
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ture. Stirring the yeast in the reaction cup to keep it well aerated was 
performed by means of a tiny magnetic flea made of a piece of fine steel 
wire sealed in a glass capillary. 

The inhibiting effects of six different wave-lengths were studied: 2345, 
2654, 2804, 3025 and 3130 A. For comparison of efficiency, determinations 
were made of the minimum doses required to prevent adaptation for four 
hours after addition of galactose. Normal adaptation time of this strain 
of yeast is 90 minutes.’ The criterion of adaptation was a measurable in- g 
crease in the rate of respiration above the endogenous before 240 minutes 
had elapsed. In doubtful cases the experiment was followed for another 
hour. Table | shows the results which were obtained, each being the re- 
sult of at least three series of experiments. 


ABSORP TION ACTION SPECTRUM 
BY NUCLEIC ACID OF ADAPTATION 
INHIBITION 


ABSORPTION 
BY 
SERUM ALBUMIN 
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RELATIVE 


WAVELENGTH IN ANGSTROMS 
FIGURE 1 


Action spectrum for inhibition by ultra-violet light of adaptation of yeast to galactose. 


Analysis of the above data reveals that there is maximum efficiency in 
the region of 2600 A. with a decrease on either side. If the relative effi- 
ciency is plotted against wave-length so that the maximum of the action 
spectrum is at the same height as that of the absorption spectrum of nucleic 
aeid, it is seen that there is excellent agreement between the two. 

Discussion.-We have seen in Fig. 1 that the absorption spectra of un- 
conjugated proteins such as serum albumin differ markedly from that [- 
of nucleic acid. Unconjugated proteins, due to the aromatic amino acids | 
they contain, have an absorption peak in the region of 2800 A.; nucleic : 
acids absorb maximally at 2600 A. because of their purine and pyrimidine 
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rings. On the basis of present evidence and theory, it would have been 
reasonable to expect that we would obtaiy an action spectrum of either, 
since both have been implicated in protein synthesis. The subject of 
synthesis of proteins, including adaptive enzymes has been reviewed by 
Northrop. He points out that enzymes, themselves proteins, must cata- 
lyze the formation of other proteins. In addition, Northrop advances the 
hypothesis that protein synthesis occurs in two steps. The first is the 
auto-catalytic formation of proteinogens or “type” proteins. The second 
step is the production of stable proteins by either catalytic or non-catalytic 
reactions. The first requires energy, the second does not. The protein- 
ogen, though hypothetical, is thought to contain nucleic acid because of 
its autocatalytic formation. Spiegelman also assigns a prominent role to 
nucleoprotein in his plasmagene theory of adaptive enzyme formation.’ 
Spiegelman and Kamen correlated the phosphate turnover of the fraction 
of yeast containing nucleic acid with adaptation.*. They stated later 
however, that their results might apply equally well to constituents in the 
fraction.” 

Although it is difficult at this time to say where nucleic acid acts in the 
production of galactozymase, it is evident from the results presented in this 
paper that it is of considerable importance. If we allow ourselves to gen- 
eralize, it would seem that all other protein synthesis would be adversely 
affected by ultra-violet light. This is not hard to imagine, considering the 
lethal effect of ultra-violet radiations; death is the logical fate of an organ- 
ism which cannot replace its outworn protoplasm. It is significant in this 
regard that the action spectrum for suppression of division of yeast and 
bacteria also resembles absorption by nucleoprotein.®  ™ 

Summary.—The action spectrum of ultra-violet light inhibition of adap- 
tation of yeast to galactose was determined. The results strongly indi- 
cate that nucleic acid is the cellular constituent affected by the light. It 
follows that nucleic acid is at least partially responsible for the adaptive 
enzyme production. The possible significance of these findings is briefly 


discussed. 


* Part of a thesis submitted as a partial requirement for the Ph.D. The author is 
indebted to Dr. Arthur C. Giese who suggested the problem and helped in its execution. 

t Present address; Dept. of Zoology, University of Massachusetts, Amherst. 
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THE STRUCTURE OF THE HETEROCHROMATIC PART OF THE 
¥-CHROMOSOME IN DROSOPHILA BUSCKI* 


By J. D. KrIvsHENKO 
DePARTMENT OF ZooLoGy, UNiverstry or Mrssour:, Mo. 
Communicated by Th. Dobzhansky, October 30, 1950 


In the salivary glands of D. busckt the Y-chromosome was at first shown 
to be a small euchromatic element whose length is about '/s of the X- 
chromosome, and consisting of some 14-15 bands."* In ganglial meta- 
phases, this Y-chromosome has the form of a tapering rod, narrow at the 
proximal end which carries a small knob-like thickening; the length is 
about equal to the Y-chromosome. This non-correspondence between 
the appearance of the Y-chromosome in the two types of cells was at- 
tributed to a heterochromatic or inert part, which is greatly reduced in 
relation to length, or apparently absent, in salivary glands. Until recently 
few attempts were made to identify such heterochromatic elements in 
salivary gland cells. 

It has been shown previously that the }-chromosome in D. buscki is 
necessary for the life of males: NXO-males die at an early stage of the 
development.*~* Asa result of a comprehensive cytogenetic investigation 
on the Y-chromosome in D. buscki, it has been possible to clarify the role 
of constituent parts in some detail, and to show their action on the develop- 
ment of males. In particular, it has been shown that the Y-chromosome is 
two-armed: the right, short arm, is the knob-like thickening at the proximal 
end, and is the euchromatic element in the salivaries; the left, long arm, is 
the inert heterochromatic part (unpublished). 

New information concerning the structure of this inert part has also been 
obtained in the course of this work, and is reported in the present paper. 

Experimental Data.—The material for this work was the mutant line 
Curled- ¥(404) obtained in an F, from crossing X-rayed males to wild type 
females. The ends of the wings of the mutant flies are bent upwards. 
This mutant behaves genetically as Y-linked, that is, it is inherited from 
father to son. Among the offspring of mutant males, about 10-12% of 
non-mutant males appear. This last class of males have a slightly short- 
ened body, and are completely sterile. 

Cytological analysis of the mutant males shows a break in region 7 of 
IIR. The proximal part of the broken IIR usually pairs fully with its 
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unbroken partner, while the distal part usually pairs only partially, or lies 
entirely free in the nucleus near the proximal ends of the other chromo 
somes. In the class of stenle males mentioned above (second class) a 
trivalent for the distal part of IIRisseen. This third, short element, which 
is the distal part of IIR seen in mutant males, sometimes pairs with its 
normal partners, or sometimes is free. The euchromatic Y-chromosome 
element in either type of male is not involved in the abberation. In 
ganghal metaphases of mutant males there is one autosome with a markedly 
longer arm and a Y-chromosome which is short and two-armed. The 
males of the second class have normal autosomes, but also a short two- 
armed Y-chromosome (these males are hy- 
poploid for a part of the Y and hyperploid 
for the distal part of IIR). 

From all these data, it is apparent that 
we are here dealing with a reciprocal trans- 
location between the Yand IIR. The posi- 
tion of the break in IIR is known; in the 
}-chromosome, the break is in the proximal 
part, near its knob-lke end. One arm of the 
new }-chromosome is this knob-like thicken- 
ing, the other arm is the distal part of the 
translocation IIR. The long heterochro- 
matic arm of the }-chromosome is translo- 
cated to the top of the proximal part of ITR, 
and the lengthened arm thus formed can be 
seen in an autosome in metaphase prepara- 
tions. This last circumstance has given us 
the opportunity to observe in salivary gland 
preparation the translocated heterochro- 
matic partof the Y-chromosome. It can be 
seen at the end of the proximal part of IIR. 
If therefore, this heterochromatic part 1s rep- 
resented by some sort of structure, this strue- 
ture must be located at the end of the proximal part of IIR. 

At first, all attempts to find any such structure were unsuccessful. 
However, trials of culture media were made to secure better development 
of D. buscki larvae, and this was responsible for final success. In one 
preparation, a unique type of chromosome element was found at the 
proximal end of IIR. The nature of this element could not be in doubt: 
it was in fact the heterochromatic part of the Y-chromosome, hitherto 
considered invisible. The most striking feature of this structure was its 
much greater transparency (as compared with the euchromatic portions), 
with pale bands, which are frequently dotted. In later work, it was often 
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Translocation involving the 
left arm of the } and the right 
arm (IER) of the second chromo 
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possible to obtain preparations showing this element, but usually the 
structure hardly stains at all and could only be seen because of differences 
in refractive index. However, among a great many preparations, a few 
were obtained which showed the entire element and many of the features 
of its structure (see figure 1). 

The length of the translocated Y-chromosome portion is about */, or 
*/, of the length of the X-chromosome, or of any of the other elements of 
the salivary gland chromosomes. Its width is about equal to that of 
ordinary paired chromosomes. The distal end, or the tip, has no re- 
semblance to any of the other salivary chromosomes of D. buscki. The 
bands are pale, difficult to see, often dotted, and all of them are single 
(dotting is not a constant character; in some preparations bands appear 
continuous). The spacing of the bands along the length of the chromo- 
some is about the same as that of the bands of the euchromatic parts of 
other chromosomes. Because of the generally uniform staining of the 
individual bands, one gets an impression that they are also more uniformly 
distributed along the length. It should be noted that observation of the 
entire length of this chromosome is very difficult, because its visibility de- 
pends on the background on which it happens to lie. The best background 
is a uniform layer of cytoplasm; in places where there is no cytoplasm, 
the chromosome cannot be seen. By the least pressure on the cover slip, 
this chromosome is very easily deformed and disintegrated, leaving no 
visible trace, indicating an exceedingly fragile structure. Presumably 
for these reasons it is extremely difficult to observe the structure of the 
Y-chromosome in individuals without the translocation. 

Discussion.—The original concept of heterochromatic chromosomes or 
chromosome parts as having no structure, concentrated as an amorphous 
mass in the so-called chromocenter, was radically modified by later investi- 
gators. In many chromosomes, genetically inert regions were demon- 
strated near the centromeres. Prokofyeva* * showed that such inert 
regions occurred not only near the centromeres but were also intercalated 
among euchromatic regions. The same author was first able to observe 
the ¥-chromosome in salivary gland cells.* In D. melanogaster it can be 
seen as a short element containing some 8-10 very thin bands. This small 
Y-chromosome is usually localized at the periphery of the chromocenter, 
and often pairs with the inert region of the centromere end of the X- 
chromosome.’ And yet, according to Prokofyeva, the observed segment 
of the Y is only its genetically active part, the remainder being invisible. 
Bauer” observed the }¥-chromosome of D. pseudoobscura as a bundle of 
chromonemata, with chromioles (the so-called loose heterochromatin). 
He was unable to see a Y-chromosome in species of Drosophila other than 
D. melanogaster and D. pseudoobscura. In D. miranda, MacKnight'! 
demonstrated the Y-chromosome as a complicated structure localized at 
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the center of the nucleus, together with the chromocenter. The pecu- 
liarity of the structure of the Y of this species is that it includes numerous 
sections of euchromatic material, originally parts of the third chromosome 
which has become involved in the sex-determining mechanism. These 
have in the course of time been much modified by fragmentation. Between 
the euchromatic segments are heterochromatic masses having no definite 
structure. In the D. nebulosa, according to Pavan,"* the V-chromosome 
is represented by a few chromomeres. 

The inert part of the Y-chromosome of D. busckt which has the definite 
structure described above is unique. First of all, the translocation having 
*/, or */, of the length of the X-chromosome makes us suppose that the 
entire }-chromosome in salivaries is proportionally as long as it is at 
metaphase; that is, the ¥ is about equal in length to the XY. This con- 
clusion conflicts with the suppositions of Muller and Gershenson'* and 
Prokofyeva’ who suggested that the metaphase Y-chromosome consists of 
heterochromatic blocks invisible in cells of the salivary glands, or alter- 
natively that it is unspiralized in metaphases. On the other hand the 
demonstration of a definite structure in the inert Y-chromosome confirms 
the conclusion of Muller and Prokofyeva' concerning the degree of re- 
semblance between euchromatic and heterochromatic segments. We 
might say that the external appearance of the heterochromatic part of the 
Y-chromosome in D. buscki seems to be a replica of a normal euchromatic 
chromosome. But these authors made conclusions about the structure 
of heterochromatic chromosomes in general, on the basis of such small 
segments as had not yet fully lost genetic activity, and which perhaps for 
that reason had not yet lost all structural organization. Similarly we 
might suppose that the }-chromosome of D. buscki is an example of such 
an element, at an intermediate stage in the process of inactivation. 

Two facts are relevant to the above problem. First, the presence in 
XX Y females of an entire }-chromosome, including its active euchromatic 
part, in no way modifies their appearance or viability. Second, males 
lacking almost the entire long left arm of the Y-chromosome are fully 
viable. Hypoploid males of our mutant line Cur/ed-Y thus behave, as 
far as the viability and the phenotype are concerned, exactly like XO, 
XY‘ and XY" of D. melanogaster. One may conclude then that the degree 
of genetic inactivation of the left arm of the Y-chromosome in D. buscki 
approaches that found in D. melanogaster, But the question still remains 
open why in D. buscki the Y-chromosome should still be visible in salivary 
gland cells even though difficult to observe, while in D. melanogaster only 
an insignificant element of it remains, and in other species are even com- 
pletely absent. 

Summary.--A translocation between the F-chromosome and an auto- 
“Gee some is described in D. busckt. The translocation provides a means of 


| 
We 
L 
| 
| 
i 
if 
; 
| 
| 
pars 
‘a \ 
1 
> 
I 
4 


Vor. 36, 1950 GENETICS: J. D. KRIVSHENKO 707 


observing in salivary gland preparations the structure of a considerable 
part of the heterochromatic region of the Y-chromosome. This inert 
segment resembles a euchromatic chromosome in size and in having clear- 
cut bands, but it is extremely transparent and takes up almost no stain. 
The bearing of this on the genetic nature of inert chromosomes is discussed. 
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EXPERIMENTS WITH THE CHEMOSTAT ON SPONTANEOUS 
MUTATIONS OF BACTERIA 
By AARON Novick AND Leo SZILARD 
Inerirure or RapromroLocy AND Broprysies, University or Carcaco 
Communicated by H. J. Muller, October 18, 1950 
Introduction.-All bacteria require for growth the presence of certain 
inorganic chemical components in the nutrient, such as potassium, phos- 


phorus, sulphur, etc., and with a few exceptions all bacteria require an 
energy-yielding carbon source, such as, for instance, glucose or lactate, 


etc. In addition to these elements or simple compounds, certain bacteria 
require more complex compounds, for instance an amino acid, which they 
are not capable of synthesizing. For the purposes of this presentation, 


any of the chemical compounds which a given strain of bacteria requires 
for its growth will be ¢alled a ‘growth factor.” 

In general, the growth rate of a bacterial strain may be within very 
wide limits independent of the concentration of a given growth factor; 
but since at zero concentration the growth rate is zero, there must of 
necessity exist, at sufficiently low concentrations of the growth factor, a 
region in which the growth rate falls with falling concentration of the 
growth factor. It therefore should be possible to maintain a bacterial 
population over an indefinite period of time growing at a rate considerably 
lower than normal simply by maintaining the concentration of one growth 
factor the controlling growth factor—at a sufficiently low value, while the 
concentrations of all other growth factors may at the same time be main- 
tained at high values. 

We shall describe further below a device for maintaining in this manner, 
over a long period of time, a bacterial population in the growth phase at a 
reduced growth rate and shall refer to it as the Chemostat. 

If the growth rate of a bacterial population is reduced, it is not a priori 
clear whether the growth rate of the individual cells which constitute 
the population is uniformly reduced or whether a fraction of the total 
cell population has ceased to grow and is in a sort of lag phase, while the 
rest keeps growing at an undiminished rate. We believe that under the 
conditions of our experiments, to be described below, we are dealing with 
the slowing of the growth rate of the individual cells rather than the cessa- 
tion of growth of a fraction of the population. 

By using an amino acid as the controlling growth factor we were able to 
force protein synthesis in the bacterial population to proceed at a rate ten 
times slower than at high concentrations of that amino acid. It appears 
that we are dealing here with a hitherto unexplored ‘‘state” of a bacterial 
population—a state of reduced growth rate under the control of a suitably 
chosen growth factor 
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The study of this “slow-growth-phase” in the Chemostat promises to 
yield information of some value on metabolism, regulatory processes, 
adaptations and mutations of micro-organisms; the present paper, how- 
ever, is concerned only with the study of spontaneous mutations in bacteria. 

There is a well-known mutant of the B strain of coli, B/1, which is re- 
sistant to the bacterial virus T,, sensitive to the bacterial virus T;, and 
which requires tryptophane as a growth factor. We used this strain and 
mutants derived from it in all of our experiments here reported. As a 
nutrient medium we used a simple synthetic lactate medium (Friedlein 
3 medium) with tryptophane added. As the controlling growth factor, we 


used either lactate or tryptophane, x 
Experiments on Growth Rates at Low Tryptophane Concentrations.—In ¢| 

order to determine the growth rate of B/1 as a function of the tryptophane q 
concentration (at high lac- } 
tate concentrations) we | 
made a series ofexperiments 
in which we incubated at 
37° at different initial tryp- 
tophane concentrations ¢, 4 | > 7 
flasks inoculated with about | 

: 100 bacteria per ce. and ob- ie | F 
tained growth curves by de- 
termining (by means of col- ad 
ony counts) the number of 
viable bacteria as a function 

of time. Because the bac- 

teria take up tryptophane, 

: the tryptophane concentra- Experiment of September 18, 1950, at 37°C. 
tion ¢ decreases during the The curve marked SLOW relates to strain B/1 

and the curve marked FAST relates to B/1/f. 


growth of the culture and 
the growth rate for the ini- 
tial tryptophane concentration c must therefore be taken from the early part 
of the growth curve. 

The growth rate a is defined by 


1 dn 
n dt 


where » is the number of bacteria per cc. 


The reciprocal value, r = ~, 
a 


we shall designate as the “generation time.’ From the generation time 
thus defined, we obtain the time between two successive cell divisions by 
multiplying by In 2. 

In figure | the curve marked “‘slow’’ shows the growth rate a as a function 
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of the tryptophane concentration ¢ for 37°. At low tryptophane concen- 
trations ¢, the growth rate at first rises proportionately with the concen- 
tration; with increasing concentrations, however, the growth rate ap- 
proaches a limit and for concentrations above 10 7/1. (micrograms per 
liter) the growth rate is no longer appreciably different from its highest 
attainable value. This highest value corresponds to a generation time of 
r = 70 min, One half of the highest value is reached at a tryptophane 
concentration of about ¢ * | y/l This concentration corresponds to 
about three molecules of tryptophane per 107" ec. (The volume of one 
bacterium is about 10~" ec.) 

The proportionality of the growth rate with the concentration of trypto- 
phane at low concentrations becomes understandable if we assume that the 
uptake and utilization of tryptophane by the bacterium requires that a 
tryptophane molecule interact with a molecule of a certain enzyme 
contained in the bacterium and that the uptake of tryptophane by these 
enzyme molecules in the bacterium becomes the rate-limiting factor for the 
growth of the bacterium. On the basis of this argument, we believe that 
down to as low concentrations of tryptophane as the proportionality of 
growth rate to concentration can be experimentally demonstrated, the 
observed growth rate of the bacterial culture represents the growth rate 
of the individual bacterium and that no appreciable fraction of the popula- 
tion goes into lag. 

The Theory of the Chemostat.--In the Chemostat, we have a vessel (which 
we shall call the growth tube) containing V cc. of a suspension of bacteria. 
A steady stream of the nutrient liquid flows from a storage tank at the rate 
of w ee.’ sec. into this growth tube. The content of the growth tube is 
stirred by bubbling air through it, and the bacteria are kept homogeneously 
dispersed throughout the growth tube at all times. An overflow sets the 
level of the liquid in the growth tube, and through that overflow the 
bacterial suspension will leave the growth tube at the same rate at which 
fresh nutrient enters the growth tube from the storage tank. 

After a certain time of such operation, at a fixed temperature, a stationary 
state is reached in the growth tube. We are interested in this stationary 
state in the particular case in which the growth rate of the bacteria is 
determined by the concentration in the growth tube of a single growth 
factor (in our specific case tryptophane). By this we mean that the 
concentration of a single growth factor (tryptophane) in the growth tube 
is so low that a small change in it appreciably affects the growth rate of the 
bacteria, and at the same time the concentration of all other growth factors 
in the growth tube is so high that a small change in them has no appreciable 
effect on the growth rate of the bacteria. As we shall show, under these 
conditions the concentration ¢ of the growth factor in the growth tube 
in the stationary state, for a fixed flow rate w, will be independent of the 
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concentration @ of this growth factor in the nutrient liquid in the storage 
tank. 

In order to see this, we have to consider the following: 

1. For zero flow rate of the nutrient (w = 0), the bacterial concentra- 


ldn 
tion m would rise in the growth tube eccording to oa = a(c), where a 


is the growth rate which, according to our ates is a function of the 
concentration, c, of the growth factor. 

2. In the absence of growth, the bacterial concentration in the growth 
tube would decrease for a given flow rate w according to the formula 


of the growth tube, 


where 3 = 8 may be called the ‘“washing-out rate” 


and 2 the washing-out time. 


After a while, for any given flow rate w, a stationary state will be reached 
in the Chemostat at which the growth rate a will be equal to the washing- 


out rate 8 (ana the generation time r equal to the washing-out time } 


i.€., 


we 1 V 
y’ 7; = 8 as (1) 
Thus, in the stationary state for any fixed flow rate w, the growth rate 
a is fixed; since a is a function of the concentration c in the growth tube, 
it follows that c is also fixed and independent of the concentration a of the 
growth factor in the storage tank. 
It may be asked what is the mechanism by which, for different values of 
a but the same flow rate w, the same concentration ¢ establishes itself in 
the growth tube in the stationary state. Clearly what happens is this: 
Suppose that, for a certain concentration a, of the growth factor in the 
storage tank, a stationary state with the concentration c in the growth tube , 
has established itself and subsequently the concentration of the growth fac- 
tor in the storage tank is suddenly raised to a higher value a;. When this 
change is made, the concentration c in the growth tube will at first rise and 
along with it will rise a, the growth rate of the bacteria, which is a function 
of c. The concentration of the bacteria in the growth tube will thus start 
to increase, and therefore the bacteria will take up the growth factor in 
the growth tube at an increased rate. As the increase of the bacterial 
concentration continues, the growth rate of the bacteria will, after a while, 
begin to fall and will continue to fall until a new stationary state is reached 


alc) = B= 
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at which the bacteria again grow at the same rate at which they are washed 
out, i.e., for which again we have a = y When this state is reached, the 


concentration of the growth factor in the growth tube is again down to the 
same value ¢ which it had before the concentration of the growth factor 
in the storage tank was raised from @, to a, while the bacterial density is 
now higher. 

In the stationary state the tryptophane balance requires that the 
following equation hold: 


azc++n Fic) (2) 
w 


azc+n (3) 
alc) 
where F(c) gives in grams per second the amount of the growth factor 
which one bacterium takes up per second. 
As can be easily seen, the amount Q of the growth factor that is taken up 
per bacterium produced is given by 
Fie) 


alc) 
so that, for the stationary state, we may also write 


awzsc+nQ or 


and for the ¢ a we may write 


(5) 

The Use of Tryptophane as the Controlling Growth Factor-—Since in the 
stationary state the tryptophane concentration in the growth tube of the 
Chemostat is always below 10 7/1. whenever the generation time is ap- 
preciably above 70 min., we may use the approximation given in equation 
(5) whenever the tryptophane concentration a in the storage tank ts above 
100 y 1. 

In order to determine the amount of tryptophane, Q, taken up per bac- 
terium produced, we grew bacterial cultures in lactate medium with varied 
amounts of tryptophane added. We found that if the initial tryptophane 
concentration is kept below 10 y/1,, then the amount of tryptophane taken 
up per bacterium produced is not dependent on the tryptophane concen- 
tration and has a value of Q = 2 X 10°'' gm. At higher tryptophane con- 
centrations, however, more tryptophane is used up per bacterium produced. 
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From equation (5), using the value of Q = 2 X 10~" gm. we obtain n = 
5 X 10’/ec. for a = 100 y/1. and we obtain » = 5 X 10°/cc. fora = 
1000 y/1. 

From this, it may be seen that, by choosing suitable values for 4 and w, 
we may vary over a wide range, independently of each other, the bacterial 
concentration » and the tryptophane concentration ¢. 

When we grew B/1 in a Chemostat (V = 20 cc.) for ten days at 37° ata 
generation time of r = 2 hrs. and at a bacterial density of 5X 10° /ce., we 
found that a change from the original bacterial strain, B ‘1, had taken place. 
The new strain, which we shall designate as B, 1/f, differs from the original 
strain only inasmuch as it grows, at very low tryptophane concentrations, 
about five times as fast as the original strain. The growth rate at higher 
tryptophane concentrations is not perceptibly different, nor could we detect 
any other difference between the two strains. The curve marked “‘fast”’ 
in figure 1 gives the growth rate of the B 1/f strain as a function of the 
tryptophane concentration at 37°. 

The ability of the B/1/f strain to grow faster at very low tryptophane 
concentrations gives it an advantage over the B/'1 strain under the condi- 
tions prevailing in the growth tube of the Chemostat; and a mutant of 
this sort must, in time, displace the original strain of B/1. 

Because in our experiments we would want to avoid—as much as 
possible—population changes of this tvpe in the Chemostat, we used in all 
of our experiments reported below this new strain, B/1/f. 

Spontaneous Mutations in the Chemostat.—If we keep a strain of bacteria 
growing in the Chemostat and through spontaneous mutations another 
bacterial strain is generated from it, then the bacterial density »* of the 
mutant strain should (for n* < m) increase linearly with time, provided 
that, under the conditions prevailing in the Chemostat, the new strain has 
the same growth rate as the original strain, so that there is no selection 
either for or against the mutant. In the absence of selection we have 


(6) 


where n* is the density of the mutant population, » is the density of the 

population of the parent strain and \ the number of mutations produced 

per generation per bacterium. Equation (6) holds under the assumption 

that back mutations can be neglected. From (6), we obtain for n* <n 
n* 


= + Const. (7) 
n r 


From this it may be seen that-~as stated above—the relative abundance 
of the mutants must increase linearly with time if there is no selection for 
or against the mutant. 
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If the growth rate of the mutant strain is smaller than the growth rate 
of the parent strain (a* < a) so that there is selection against the mutant 
in the growth tube of the Chemostat, then the density n* of the mutant 
population should—after am initial rise—remain constant at the level 
given by 


(8) 


Experiments on Spontaneous Mutations in the Chemostat,—Of the various 
mutations occurring in a growing bacterial population, mutants resistant 


ff 


FIGURE 2 FIGURE 3 


Experiments of May 3, 8, and 28, Experiment of July 19, 1950, at 37°C 
1950, at 37°C. giving for strain giving for strain B/1/f the concentration 
B/1/f for three different values of of mutants resistant to Ts; (left-hand 
the generation time the concentration scale) and mutants resistant to T, 
of the mutants resistant to T,, for (right-hand seale) for a population 
a population density of 5 x 10° density of 25 K 10° bacteria perce. In 
bacteria per cc this experiment oxygen containing 0.25% 

CO, was used for aeration 


to a bacterial virus are perhaps the most easily scored with considerable 
accuracy. In our experiments we mostly worked with mutants of our 
coli strain which were resistant to the bacterial viruses Ts or Ty. 

When we grow the strain B/1/f in the Chemostat with a high concen- 
tration of tryptophane but a low concentration of lactate in the nutrient 
in the storage tank, so that lactate rather than tryptophane is the con- 
trolling growth factor, we find-—after a short initial period—that the 
bacterial densities of the mutants resistant to Ts or Ts each remain at a 
constant level. These levels appear to correspond to a selection factor 
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~~ of a few per cent. 


We are inclined tentatively to assume that the behavior of these two 
mutants exemplifies the general rule that the vast majority of all the 
different mutational steps leading away from the wild type yield mutants 
which—under conditions of starvation for the carbon source—-grow slower 
than the parent type. 

On the other hand, if we grow our tryptophane-requiring strain in the 
Chemostat with a high concentration of lactate but a low concentration 
of tryptophane in the nutrient in the storage tank (so that tryptophane 
rather than lactate is the controlling growth factor) and if we run the 
Chemostat at a generation time well above 70 min. (the generation time 
at high tryptophane concentrations)—-then there is no reason to expect 
mutants in general to grow appreciably slower than the parent strain, 
particularly if the growth of the parent strain is kept very slow by keeping 
the tryptophane concentration in the growth tube very low. In this 
case one would rather expect a mutation to affect the growth rate only if 
it affects the uptake or utilization of tryptophane by the bacterium or if 
the mutant is a very slow grower. Accordingly, we should, in general, 
expect the mutant population to increase linearly with time in the Chemo- 
stat when tryptophane is used as the controlling growth factor. 

Figure 2 gives for 37° the experimental values for the bacterial density for 
the mutant population resistant to T, in the growth tube of the Chemostat 
as a function of the number of generations through which the parent strain 


t 
has passed in the Chemostat. (Nuniber of generations g = ‘) The 
T 


three curves in the figure correspond to generation times of 2 hours, 6 hours 
and 12 hours. The slope of the straight lines gives A, the mutation rate 
per generation, as 2.5 X 10~*; 7.5 KX 10-*; and 15 X 10~* per bacterium, 
We see that the mutation rate per generation for r = 6 hours is three times 
as high and for r = 12 hours is six times as high as it is for r = 2 hours. 
Thus the mutation rate per generation is, in our experiment, not constant 
but increases proportionately with r and what remains constant is the 
number of mutations produced per unit time per bacterium. According 


to the above figures, we have ~ = 1.25 X 10~* per hour per bacterium. 
Tt 


This result is not one that could have been foreseen. If mutants arose, 
for instance, as the result of some error in the process of gene duplication, 
then one would hardly expect the probability of a mutation occurring per 
cell division to be inversely proportionate to the rate of growth. 

If the processes of mutation could be considered as a monomolecular 
reaction—as had been once suggested by Delbruck and Timofeeff- 
Ressovsky~—then, of course, the rate of mutation per unit time should be 
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constant. The rate & of a monomolecular reaction is given by 


kh = (9) 


The value of the constant A can be calculated from the observed reaction 
rate & and the heat of activation W (which can be obtained by determining 
the temperature coefficient of the reaction). 

Using the Chemostat, we have determined the rate of mutation to 
resistance to T, at 25° (for r = 6 hrs. and r = 12 hrs.) and found it to be 
about one half of the mutation rate at 37°. From this value and the 


mutation rate of = 125 % 10°°% per hour per bacterium at 37° we 


compute A per sec. 
In a condensed system, such as an aqueous solution, A has been found 


to lie between 10° and 10'* per see. for known monomolecular reactions. 


Therefore if the mutation studied by us were due to a monomolecular 
reaction, it would have an A value 10° times lower than the lowest value 
so far found, 

The density of the mutants resistant to the bacterial virus Ts, in the 
Chemostat, with tryptophane as the controlling growth factor, also appears 


to rise linearly with time for r = 2 hours, r = 6 hours and r = 12 hours, 


but our results so far are not sufficiently accurate to say whether this 


mutation also occurs at a constant rate per unit time for different genera- 
tion times r. The temperature coefficient of the mutation rate appears 
to be very low, but again this conclusion must await more accurate experi- 


ments 

Phe result obtained for mutation to resistance to the virus Ts, showing 
that this mutation occurs at a constant rate per unit time up to a generation 
time of r = 12 hours, raises the question whether this is generally true of 
spontaneous bacterial mutations or whether we are dealing in our case 
with certain exceptional circumstances. Clearly, a number of different 
mutations will have to be examined, different amino acids will have to be 


used as the controlling growth factor and other conditions will have to be 
varied before one would draw the far-reaching conclusion that our observa- 
tion on mutation to resistance to the virus T, exemplifies a general rule. 
Vutants Resistant to T,.—We find that mutants resistant to Ty are 
selected against in the Chemostat when grown either with lactate or with 
tryptophane as the controlling growth factor, i.e., the number of mutants 


remains-—after an initial rise—at a fixed level. 

It is known that of the different mutants of the B strain of coli which are 
resistant to the virus Ty, the most frequent one is also resistant to the 
viruses T, and T; and that this mutant is a very slow grower under ordinary 
conditions of culture. It is conceivable that this might explain why the 
mutants resistant to T, are selected against in the Chemostat even when 
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the bacterial population grows under tryptophane control and at a much 
reduced rate. 

Manifestation of ‘Evolution’ in the Chemostat.—lf a bacterial strain is 
grown over a long period of time in the Chemostat, from time to time a 
mutant might arise which grows faster, under the conditions prevailing 
in the Chemostat, than the parent strain. If this happens, practically 
the entire bacterial population in the Chemostat will change over from 
the parent strain to the new strain. We have discussed one change-over 
of this sort, i.e., the change-over from the strain B/1 to the strain B/1/f. 
There is no reason to believe, however, that no further change-over may 
take place when we start out with B/1/f as the parent strain and continue 
to grow it in the Chemostat over a long period of time, 

We have seen that the mutants resistant to T; accumulate in the Chemo- 
stat and that their number rises linearly with the number of generations, 
giving a straight line, the slope of which is given by A. If now at a certain 
time the population changes over in the Chemostat from the parent strain 
to a faster-growing strain, the accumulated mutants resistant to the 
bacterial virus T; which were derived from the parent strain should dis- 
appear from the Chemostat along with the parent strain. This should 
lead to a fall in the number of mutants resistant to the bacterial virus T, 
during a change-over from the parent strain to the faster-growing strain. 
After the change-over to the new strain, the concentration of the mutants 
resistant to T, may be expected again to increase linearly with the number 
of generations, giving a straight line which has the same slope as before 
the change-over, because the new strain which displaces the parent strain 
may be expected to mutate to resistance to T, at an unchanged rate A. 

Thus, we may in general expect, when a change-over in the population 
takes place, the covcentration of the mutants resistant to Ts to shift from 
one straight line which lies higher to another, which lies lower. The 
magnitude of this shift may be somewhat different from experiment to 
experiment, depending on when mutants resistant to Ty happen to make 
their first appearance in the population of the new strain. 

At the outset, the bacteria belonging to the new strain will be few in 
number but their number will increase exponentially with the number of 
generations until-——at the time of the change-over——the bacteria belonging 
to the new strain become an appreciable fraction of the total population. 
If the mutation rate to resistance to T, is of the order of magnitude of 
10°, then it is unlikely that such a mutant should appear in the population 
of the new strain until its population has reached perhaps 10’. However, 
because an element of chance is involved, occasionally a mutant resistant 
to T; may appear earlier and, if that happens, the “‘shift’’ associated with 
the change-over will be smaller and in principle it might even be negative. 

If a bacterial population remains growing in the Chemostat for a suffi- 
ciently long time, a number of such change-overs might take place. Each 
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such step in the evolution of the bacterial strain in the Chemostat may be 
expected to manifest itself in a shift in the ascending straight line curve of 
the T,; resistant mutants. 

As we have seen, the mutants resistant to T, remain—apart from an 
initial rise--at a constant level in the Chemostat. However, when the 
bacterial population in the Chemostat changes over from a parent strain 
to a new strain, the T, resistant mutants might change over from one level 
to another, because the selection against the two strains might be different. 

Figure 3 shows, for mutants resistant to T, and for mutants resistant to 


t 
T,, the number of mutants as a function of the number of generations - 
T 


in a Chemostat which was run for 300 hours at r = 4 hours with trypto- 
phane as the controlling growth factor. 

It may be seen that these two curves show a population change-over of 
the type just described. The curve for the T; resistant mutants shows a 
shift, P, of P = 32 generations. 

A number of shifts of this type were observed in different experiments. 
We verified that these “‘shifts’ represent population change-overs by 
showing in one case that (under the conditions prevailing in the chemostat) 
bacteria taken from the Chemostat before the change-over in fact grow 
slower than bacteria taken from the Chemostat after the change-over. 

In order to show this, we took from the Chemostat before the change- 
over a bacterium resistant to T; and after the change-over a bacterium 
sensitive to T, and inoculated a second Chemostat (operated under identical 
conditions) with a 50-50 mixture of these two strains. We then found that 
the relative abundance of the resistant strain rapidly diminished. In the 
corresponding control experiment we took a sensitive bacterium from the 
Chemostat before the population change-over and a resistant one after 
the population change-over and again found that the strain prevalent 
before the change-over (this time the sensitive one) was the slower grower. 

In the later stages of the change-over the concentration x of the original 


t 
strain falls off exponentially with the number of generations, g = -, so that 
Tr 


we may write x = Ce~"". In our experiment we obtained for y a value 
of y = 3.25. ; 

It should be noted that the value of y can be read also directly (though 
not accurately) from the curve, which gives the concentration n* of the 
resistant mutants as the function of g, the number of generations. During 
the change-over the concentration ¢ of the tryptophane in the growth 
tube goes over from an initial value ¢, to a final, lower value c, and it can 

+ 


be shown that for the midpoint of the change-over at which ¢ = 


we have 
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(10) 


where P is the magnitude of the shift expressed in the number of genera- 
tions by which the ascending straight line of the resistant mutants is 
shifted in the change-over. This formula holds only if r is large so that 
the rate of growth of the bacteria in the Chemostat is proportionate to the 
tryptophane concentration c. Because the exact position on the curve 
of the midpoint of the change-over on the curve n* is not known, this 
formula can give only a rough indication for the value of y. 

In our case, the estimate based on it gave for y a value of y = 2.4 in 
place of the directly observed value of y = 3.25. Within the limits of the 
accuracy of our curve for n* these two values are consistent with each other. 

Population change-overs manifesting themselves in a shift in the 
ascending straight line of the T, resistant mutants occurred in every ex- 
periment carried at r = 4 hrs. beyond the 50th generation: In an experi- 
ment carried to the 450th generation at a bacterial density of 2.5 X 10*/cc., 
a number of such shifts occurred, the last one at about the 350th genera- 
tion. (In the course of this experiment the mutants resistant to Ty rose 
twice from a low level to a high peak, the first of which reached 4.6 K 10* 
and the second 4.5 X 10* mutants perce. This phenomenon is now being 
investigated.) 


It may be said that our strain, if grown in the Chemostat at low trypto- 
phane concentration for a long period of time, undergoes a number of 
mutational steps, each one leading to a strain more “‘fit’’ than the previous 
one, and that each step in this process of evolution becomes manifest 
through the shifts appearing in the curve of the mutants resistant to Tg. 


THE ANALYSIS OF A CASE OF CROSS-STERILITY IN MAIZE 
By Drew SCHWARTZ 
DEPARTMENT OF Borany, UNtversity or ILLINOIS 
Commur'-ated by M. M. Rhoades, October 25, 1950 


The male gametophyte of maize is extremely sensitive to chromosomal 
unbalance. Duplications and deficiencies are almost always accompanied 
by reduced pollen transmission. The functioning of the male gametophyte 
is also known to be affected by a number of genic factors. The most 
thoroughly studied locus and probably the most interesting is Ga; on 
chromosome 4. First detected by Correns* because of aberrant F; ratios 
for sugary-starchy, this factor was subsequently studied by Jones,* ” 
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Emerson,* ’ Mangelsdorf,” and Demerec.‘ Emerson established that the 
differential fertilization was not due to the sugary gene itself but to some 
other gene with which it was linked. 

Ga, is comparable to the self-sterility genes found in strictly cross-ferti- 
lized plants such as Nicotiana, Oenothera, Trifolium, ete., in that there 
is interaction between the stylar tissue and the male gametophyte. When 
the female parent is homozygous recessive ga), ga; pollen can compete 
successfully against Ga, and fertilize half of the ovules. However, if the 
plant used as the egg parent is heterozygous or homozygous Ga;,, the ga; 
pollen is a poor competitor and achieves fertilization in only 0-4% of the 
ovules. In the absence of competition, that is, if only ga; pollen is used 
in the cross, full seed set is obtained regardless of the genotype of the female 
parent.’ Thus the incompatibility between ga; pollen and Ga, silk is not 
detectable in the absence of competing Ga, pollen. The purpose of this 
paper is to present evidence for the existence of a third allele at this locus. 


TABLE 


Summary or Data SHowinG tus Comperrrive FUNCTIONING OF Ga;", Gai, AND ga; ON 
Sry_es or Various GENoryPes 
- NUMBER OF KERNELS 

CONBTITUTIONS STARCHY SUGARY TOTAL SUGARY 
Ga;'Su/Gaysu X 19,168 5865 25,033 23.4 
Ga;'Su/gasu 4,104 1300 5,404 23 
Gasu/Gaysu X 781 1,447 46 
Gaysu/ga,Su X GaySu/Gaysu 4,624 1432 6,056 23 
Su/Gaysu Gay Su/gaysu 1,548 221 1,769 12.5 
Ga;'’Su/gajsu Ga;'Su/gaysu? 4,142 611 4,753 12 
Gaysu/Gaysu X 237 323 26 
Gaysu/ga,Su X 3,866 675 4,541 14 
ga,su/gaysw Ga, 412 408 820 49 
Su/gajsu Gaysu ga, 2,006 1104 3,200 34 


This allele, designated Ga,’ since its effect is stronger than Ga,, was found 
by Dr. M. M. Rhoades who turned it over to the author for further analysis. 

The most important difference between Ga, and Ga,’ is that ga, pollen 
completely fails to function on styles homozygous for Ga,’ even in the ab- 
sence of competing pollen. From 14 such crosses only one seed was ob 
tained; this could have been the result of contamination. Crosses of 
ga, pollen on styles heterozygous for Ga,’ yield a partial seed set. Less 
seed is set when the female parent is heterozygous Ga,’Ga, than when a 
heterozygous Ga,'ga, plant is used. The total amount of seed produced 
in such crosses varies greatly from plant to plant. 

The linkage relation of Ga,’ with Su and its interaction with Ga, suggest 
that itis an allele of Ga,. Both Ga, and Ga,’ exhibit selective fertilization in 
competing with ga, only when the female parent in the cross carries one of 
the dominant alleles. No selective fertilization occurs when the female 
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parent is homozygous ga; (table 1, item 9). On the assumption that Ga; 
and Ga,’ represent two independent loci, and that ga, and ga’ gametes can 
compete successfully on respective recessive styles against pollen carrying 
the dominant allele, no differential fertilization should result from crosses 
of ga,ga;Ga'Ga* by Gayga,ga’ga’® and of GaGayga’ga* by 
Moreover, pollen from Ga,Ga,ga'ga’ individuals should give no seed set 
when crossed on gaga;Ga'Ga’ plants. This, however, was not found to 
be true; in the crosses of Ga;"ga; X Gayga,* and Ga;'Ga,; X Ga,'ga,;*, gay 
pollen is a very poor competitor (table 1). The cross Ga;'Ga,;’ X GayGa," 
vielded a full seed set. These data negate the hypothesis that Ga,’ and 
Ga, represent two independent loci but do not rule out the possibility of an 
interaction between the two. One might propose that the presence of Ga, 
in the Gaga’ gamete permits the functioning on Ga’ styles of an otherwise 
handicapped pollen grain. It seems unlikely, however, that such inter- 
action of loci would exist in view of the specific action of the self-sterility 
alleles in Nicotiana. East and Yarnell* have shown that there are at least 15 
alleles at the self-sterility locus. Each allele is so specific in its action that 
the rate of growth of a pollen tube containing any one of them is determined 
solely by the presence of that particular allele in the diploid sporophytic 
tissue of the style and is completely independent of any other which may 
also be present. This specificity breaks down to some extent in polyploids 
of Oenothera’ and white clover.' 

Linkage studies have shown that the recombination values obtained for 
Ga,’ and Su are very close to those found by Emerson for Ga, and Su. A 
recombination value of 30% between the latter was calculated on the basis 
of 0% functioning of ga, pollen. This value was obtained by determining the 
percentage of sugary kernels resulting from a cross such as Ga,su/Gajsu X 
Ga,Su/gasu*. If no ga, pollen functions, all of the sugary kernels result 
from crossing over between Ga, and Su, and the percentage of sugary kernels 
on the ears represents the crossover frequency. However, if some ga; 
pollen does function, then not all of the sugary kernels are crossovers and 
the recombination value is somewhat lower. Emerson showed that ga; 
pollen is responsible for the fertilization of from 0 to 9% of the ovules with 
the mean at 4°). On the basis of 4° functioning of ga; pollen he calcu- 
lated a recombination value of 27.89%. When the Ga,’ allele was sub- 
stituted for Ga, in similar crosses, a percentage of 25.8% sugary kernels 
was obtained (table 1, item 6). This value, which represents the recom- 
bination between Ga," and Su, agrees well with that calculated by Emerson 
for Ga,-Su. 

That the discrepancy between the recombination values of 30% for 
Gay-Su and 25.8% for Ga;’-Su is due to the amount of ga; pollen which func- 
tions is best shown by the following experiment. Ga,'Su/gasu plants 
were used as the pollen parents in crosses with (1) Ga,su/ga,Su, and (2) 
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(ra;'Su/ga,su plants. Since the same male parent was used in both crosses 
the amount of crossing over between Ga,’ and Su would be identical and 
any significant difference in the percentage of sugary kernels resulting from 
these crosses would be due to the functioning of ga, pollen. The data are 
given as items 6 and Sin table |. The crosses involving Ga; gave 14.9% 
sugary kernels while 12.9) sugary seeds were obtained when Ga,’ was in- 
volved. The difference between 12.9 and 14.9% is significant at the 1% 
level and may logically be ascribed to the differential functioning of ga; 
pollen on Ga, and Ga,’ silks. The value of 14.9% is very close to the 
15.1%) found by Emerson in similar crosses using Ga,. It follows that ga; 
pollen occasionally functions on Ga, styles even when competing against 
pollen. 

The close agreement in the amount of crossing over between Ga," and 
Ga, with Su is indicative of an allelic relationship between the two genes. 
However, a possible alternative must be considered—namely, that Ga,’ 
and Ga, represent two loci lying equidistant from (in terms of crossover 
units), but on opposite sides of, Su. Such a condition would place the 
two genes more than 50 map units apart. 


TABLE 2 
DisrripuTion or Ears ResuLTInG THE OF PLANTS Derivep 
THE Cross Ga,'Su/Gajsu X gaysu? 


SUGARY EERNELS, % 
VAMILY 8-12 13-17 18-22 23-27 


1360 0 4 3 
1361 
1362 


In order to test this hypothesis the following crosses were made: 
(Ga,'Ga,'SuSu * Ga,Gasusu”) gaygasusu’. The sugary seeds were 
discarded since a heterozygous Su/su condition is necessary to detect 
the presence or absence of the gamete factor. On the assumption that 
two loci are involved, the F; would be genotypically ga,SuGa’/Ga,suga’. 
The second cross would yield four classes of starchy kernels: (1) ga;SuGa’ / 
g2,suga’, (2) Ga,SuGa" /ga'suga,, (3) ga,Suga’/ga’suga,, and (4) Ga,Suga*/ga’ 
suga,. The first class is a non-crossover; the second, a single crossover 
in region (1) (Ga,-Su); the third, a single crossover in region (2) (Su-Ga’); 
and the fourth represents a double crossover in regions (1) and (2). Classes 
1, 2 and 4 are difficult te distinguish and when self-pollinated would yield 
only 12-15°% sugary seeds due to the linkage of Su with the dominant 
gamete factor. However, class 4, comprising almost 25° of all the starchy 
seeds would be detected since in the absence of either Ga; or Ga’ normal 
Mendelian ratios of starchy-sugary kernels, i.e., 25% sugary, should result 
upon self-pollination. On the other hand, if Ga" and Ga,, are alleles of 
the same locus, the cross of the /; with ga,su should give only 2 classes of 
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starchy seeds, Ga;’Su/gasu and Ga,Su/gasu. On this basis, no 3:1 
starchy-sugary ratios would be expected on self-pollinated ears. The re- 
sults are summarized in table 2. Since only low sugary ratios resulted, it 
is apparent that Ga,‘ and Ge, are not separate loci located on opposite sides 
of Su. 


When Ga, and Ga,’ pollen compete, the latter has a definite advantage ' 


in achieving fertilization on both Ga, and Ga,’ styles. Tests to determine 
the degree of selective fertilization were made using female parents of four 
genetic constitutions: (1) Ga;'Su/Ga,su, (2) (3) Gaysu/Gaysu 
and (4) Gasu/ga,Su. The pollen parents were Ga,su/Ga,;'Su. The re- 
sults are given in table 1. In the absence of any differential fertilization 
the crosses with the heterozygotes should produce 25° sugary kernels and 
the backcrosses to homozygous sugary plants should yield equal numbers 
of both types of seeds. However the 23.59% sugary kernels in the F; 
population and the 46°% sugary in the backcrosses were significantly dif- 
ferent at the 1% level from the expected percentages. Assuming 26°% re- 
combination between Ga, and Su it may be calculated that Ga,’ pollen fer- 
tilized 56% of the ovules in competition with Ga, pollen. 

The manner in which selective fertilization is accomplished has not been 
established. Jones,’ working with Ga,, compared the ratios of sugary 
to starchy kernels on the butt and tip ends of the ears and concluded that 
the causal agent was differential pollen-tube growth. Emerson’ was not 
able to confirm this. Demerec* claimed that selective fertilization was due 
to cross sterility rather than to competition between the pollen tubes. 
However, from his data it would appear that Demerec was working with a 
different allele from that studied by Jones and Emerson. Demerec’s 
gamete factor is analogous to the Ga,’ allele reported in this paper in that 
ga, pollen does not function well in styles which carry the dominant allele 
even in the absence of competing pollen. The 12.4% sugary kernels which 
he obtained in the F,; of sweet corn by plants showing the cross-sterility 
compares well with the reported 12.9%) sugary resulting from the self- 
pollination of Ga,’Su/ga,su plants. 

Demerec listed three possible explanations of this cross sterility. 

(1) Inability of pollen with recessive gametophyte factor to germinate 
on silks having dominant allele. 

(2) Inability of recessive pollen tubes to grow in silks with dominant 
allele. 

(3) Inability of recessive pollen tubes to reach the ovules in the plant 
having the dominant allele. 

In the case of Ga,’, it was possible to eliminate the first two causes experi- 
mentally. Plants of the constitution Ga,;'Ga,’ were pollinated with ga, 
pollen. Two hours after pollination the silks were removed, fixed in modi- 
fied Carnoy solution and stained with propionic carmine. It was evident 
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that the pollen grains had germinated and that the pollen tubes had pene- 
trated the stylar tissue. 

Demerec’s third suggestion remains as a possibility and a fourth may 
be added which is likewise untested. The failure of ga, pollen to fertilize 
Ga,’ plants may be independent of stylar interaction altogether, and may 
involve some abnormal reaction within the ovule, such as inability of the 
pollen tube to penetrate the embryo sac. None of these explanations, 
however, are applicable in the case of Ga,-ga, competition since, in the ab- 
sence of Ga, pollen, ga, gametes function normally on Ga, styles. 

A number of cases of cross-sterility in maize may be due to the action 
of Ga,’. This is especially true of the cross-sterility found among varieties 
of popcorn.’ Although the point of origin of the Ga; allele is not known, 
the fact that both Ga; and Demerec’s cross sterility factor were found in 
popeorn supports this contention. 

' Atwood, S. S., and Brewbaker, J. L.., Ree. Genet. Soc. Am., 19, 87 (Abstract) (1950). 

* Brunson, A. M., and Smith, G. M., 7. Am. Soc. Agron., 37, 176-183 (1945). 

* Correns, C., Ber. deut. bot. Ges , 20, 159-172 (1902). 

Demerec, M., Z. 7. A. V., $0, 281-291 (1929) 

* Kast, M., and Yarnell, S. H., Genetics, 14, 455-487 (1929). 

* Emerson, R.A., Anat. Rec., 29, 136 (Abstract) (1925). 

’ Emerson, R. A., Genetics, 19, 137-156 (1935) 

* Jones, D. F., Proc. Nati. Acap. Sct., 10, 218-221 (1924). 

* Lewis, D., Biol. Rev. Cambridge Phil. Soc., 24, 472-496 (1949) 

Mangelsdorf, P. C., and Jones, D. F., Genetics, 11, 423-455 (1926). 


THE USE OF SODIUM NUCLEATE IN THE STUDY OF THE 
MUTAGENIC ACTIVITY OF ACRIFLAVINE IN ESCHERICHIA 
COLI* 

By Evertyn M. Wirxryt 
DerarrMenr or Generics, INsriruTton, Sprinc Harsor, L. N. Y. 


Communicated by M. Demerec, October 25, 1950 


In a previous study,' neutral acriflavine was shown to induce mutations 
from phage-sensitivity to phage-resistance in Escherichia coli. The pres- 
ent report concerns an improved method of investigating the mutagenic 


potency of certain compounds, using acriflavine as a model. 

One of the difficulties encountered in working with acriflavine, as well 
as with numerous other compounds, was the tenacity with vitich the 
chemical remains bound to the treated cell. Ordinary washin, methods 
were effective in removing acriflavine only to a limited degree, and thor- 
oughly washed suspensions of treated bacteria were found to contain 
enough active chemical to limit subsequent growth, and to interfere with 
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the process of infection by bacteriophage. Another problem was the 
tendency of acriflavine to cause clumping of bacteria in dense cell suspen- 
sions. The agglutination could be controlled in part by adjustment of 
experimental conditions, but was extremely difficult to eliminate entirely. 
Since these complications, which were encountered with several other 
compounds, could operate to simulate mutagenicity, it seemed important 
to develop a method whereby they could be circumvented or overcome. 

It has been reported that the antibacterial activity of acridines can be 
inhibited by the addition of nucleic acids and some of their derivatives,* * 
and the kinetics of this interaction have been described.‘ In these experi- 
ments, sodium ribonucleate was used to bring the exposure of bacteria to 
acriflavine to a sharp end point, to eliminate residual activity of this sub- 
stance in the subsequent handling of the treated bacteria, and to prevent 
agglutination. In principle, this method is applicable to other antibac- 
terial compounds for which more or less specific reversing agents are avail- 
able. 

Material and Methods.—-Strain B/r of E. coli was used. The stock was 
carried on nutrient agar slants, and subcultured every two months. Inoc- 
ula for experimental cultures were taken from a 24-hour culture grown in 
nutrient broth from about 100 cells. 

Bacteriophage Tl was used in the isolation of mutants, and was kept in 
broth suspensions having a titer of about 10" particles per ml, 

Difco nutrient broth and Difco nutrient agar, each with 0.5% NaCl 
added, were used us culture media. 

Neutral acriflavine was obtained from the National Aniline Division of 
the Allied Chemical and Dye Corporation, New York City. Solutions were 
made in distilled water, sterilized by boiling and used for no longer than 
three days. Sodium ribonucleate (yeast), obtained from Schwartz Labo- 
ratories, New York City, was’ dissolved in distilled water and sterilized 
by immersion in a boiling water-bath for fifteen minutes. 

All assays were made by quantitative dilution and surface plating. 
Determinations of the frequency of phage-resistant mutants were made 
by spreading 0.1-ml. aliquots of culture on the surface of plates previously 
coated with about 10° phage particles. 

Incubation was at 37°C. 

Each experimental series consisted of 20-50 tubes containing 5 ml. of 
broth to which acriflavine was added to give a concentration of 0.01%. 
Each tube was inoculated with about 10° bacteria and incubated for 3-4 
hours, or until. survival would be expected to reach between 107° and 

10>. At this time, sodium nucleate was added to each tube, to give a 
concentration of 0.5°), and several tubes were assayed immediately to 
determine precisely the number of viable bacteria per culture. The cul- 
tures were then incubated for 18-24 hours, after which several were assayed 
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to determine the average titer, and an aliquot from each culture was plated 
with phage to obtain the number of T1-resistant mutants. 

Controls consisted of 20-50 tubes containing 5 ml. of 0.01 acriflavine in 
broth, to each of which were added simulianeously sodium nucleate (final 
concentration 0.5%) and an inoculum of 10*-10* untreated bacteria. These 
cultures were incubated 18-24 hours, and assayed for titer and phage-resist- 
ant mutants in the same manner as the experimental cultures. 

Thus the experimental and control cultures were grown from inocula of 
comparable size in a medium containing both acriflavine and nucleate. 
In the case of the experimental cultures, the inoculum consisted of sur- 
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TIME IN HOURS 
FIGURE 1 
Killing and Growth in Experimental and Control Cultures. A-A,; = killing curve 
of B/r in 0.01% acriflavine in broth; A, = growth of acriflavine-treated bacteria after 
addition of 0.5% sodium nucleate (experimental cultures); B = growth of untreated 
bacteria in broth containing 0.01% acriflavine + 0.5% sodium nucleate (controls). 


vivors of a prior exposure to acriflavine alone, whereas in the controls un- 
treated bacteria were used. 

Experimental Results.—-Growth in Nucleate-Inactivated Acriflavine: Fig- 
ure | shows the course of killing and growth in the experimental cultures, 
and of growth in the controls. Killing proceeds logarithmically in acri- 
flavine until the addition of nucleate, at which time it is sharply arrested. 
The growth of the survivors in the nucleate-acriflavine mixture is com- 
parable to that of the controls, except for the prolongation of the lag 
phase. The final titers of experimental and control cultures were not 
significantly different. 
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Frequency Distribution of Phage-Resistant Mutants in Experimental and 
Control Cultures.—Figure 2 shows the frequency distribution of mutants 
resistant to bacteriophage T1 in 300 experimental cultures and 300 control 
cultures. The spread observed in the control series is a function of the 
spontaneous mutation rate and the clonal grouping of mutants in liquid 
cultures,’ and has been found to be extremely reproducible. It will be 
seen that the distribution of mutants in the experimental cultures is very 


columns) and 300 control cultures (white columns), 


different. The largest class of cultures in the experimental series is that 
having between 10* and 10° mutants per 10° bacteria, whereas the largest 
control class has less than 10 mutants per 10° bacteria. The frequency 
distributions for both experimental and control culture series are highly 
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MUTANTS PER 10° BACTERIA 
FIGURE 2 


Frequency distribution of Tl-resistant mutants for 300 experimenta) cultures (black 


reproducible under the conditions specified, as will be shown below in an- 
other connection. 

Selection.-The observed shift in the frequency distribution of mutants 
in experimental cultures can be accounted for by induced mutation only if 
two possible modes of selection can be ruled out. Selection could enter 
through differential survival of spontaneous mutants in the experimental 
cultures during exposure to acriflavine, or through differential stimulation 
of the rate of growth of phage-resistant mutants by the acriflavine-nucleate 
mixture, or through a combination of the two. These possibilities were 
tested in the following way: Ten single-colony isolates of phage-resistant 
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mutants were made, each colony coming from a plating on phage of a dif- 
ferent experimental culture, If selection were responsible for the increased 
frequency of mutants in these cultures, mutants thus obtained should pos- 
sess the particular advantages involved, having been through the selective 
sieve of the experimental procedure. After repeated streaking and single- 
colony isolation to ehminate contaminating phage particles, the ten phage- 
resistant strams were grown in broth. Mixtures were made of each mu- 
tant strain with the phage-sensitive parent strain, and each mixture was 
assayed to determine the relative proportions of its components. The mix- 
tures were then subjected to the experimental procedure deseribed above: 
inoculation into acriflavine, and incubation for four hours followed by addi- 
tion of sodium nucleate and incubation for 18-24 hours. Assays to deter- 


TABLE 


Seceerion Agrrercia Mixtures or 10 InpEPpENDENT T1-ResistaANT Muranr 
Sreatns B/r, Sunrecrep ro Tor EXPERIMENTAL PROCEDURE 


percent Initial proportion of Tl -resistant mutants 

VW,. per cent Proportion of Tl-resistant mutants immediately after the addition of 
wxiium sucheate 

Wy, percent Final proportion of Tl-resistant mutants 


Me Me 


19 j 6 
18.1 
1.3 


33.2 35 56 6 42.0 


mine the proportions of phage-resistant and phage-sensitive components 
were made immediately after the addition of nucleate, and ajrain after 
the I8- to 24-hour incubation period, The results are shown ta table 1. 
It is evident that selection does not favor the phage-resistant mutants, 
either by differential survival or by differential growth. 

Vutagenic Effect of Acriflavine on an Acriflavine-Resistant Strain. 
The growth of the normal B r strain is completely inhibited by a concen- 
tration of 0.002%, of neutral acriflavine in broth. By serial transfer in 
increasing concentrations of acriflavine, a strain was derived from B/r which 
could grow in 30 times this concentration, or 0.06%. This acriflavine- 
resistant strain, when treated with acriflavine according to the standard 
procedure described above, using a concentration of 0.01% acriflavine, 
gave a frequency distribution of phage-resistant mutants identical with 
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that of the untreated controls. The four-hour exposure to this concentra- 
tion of acriflavine was entirely non-toxic for the resistant strain. When the 
experiment was repeated, however, using a concentration of 0.3%, which 
gave a survival of 10~* to 10~*, the frequency distribution of phage-resist- 
ant mutants was essentially the same as that obtained with strain B/r at 
the 0.01%) concentration, which yields the same survival during the 4-hour 
treatment. Figure 3 shows the frequency distributions of pkage-resistant 
mutants obtained for the acriflavine-resistant strain, using 0.3% acri- 
flavine, as compared with those obtained for B/r at a concentration of 
It is evident that the mutagenic effect 


0.01%, and for untreated controls. 
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Frequency distribution of Tl-resistant mutants for 100 experimental cultures of 
acriflavine-resistant bacteria (cross-hatched columns); 100 experimental cultures of 
and 100 controls (white columns). 


acriflavine-sensitive bacteria (black columns) ; 


of acriflavine is related to the survival levels of the sensitive and resistant 
strains, rather than to the absolute concentration of acriflavine used in the 
treatment. 

Discussion.—The mutagenicity of acriflavine is confirmed in these ex- 
periments, under conditions free of the complications encountered earlier 
and enumerated above. The method is essentially qualitative, since the 
cultivation following the acriflavine treatment introduces clonal variance, 
which makes precise determination of the frequency of induced mutations 
extremely difficult. In addition, the method as described does not permit 
distinction between ‘‘zero-point’’ and delayed mutations, and both types 
could contribute to the greater frequency of cultures having high numbers 
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of mutants. Modifications that would permit the separation of the two 
types of induced mutations are theoretically possible, however. One 
such procedure would involve adding bacteriophage to the nucleate-inac- 
tivated acriflavine tubes mmmediately after adding nucleate, so that only 
tubes containing at least one zero-point mutant would develop turbidity. 
Addition of bacteriophage after various incubation periods would similarly 
permit the measurement of delayed mutations. Preliminary experiments 
along these limes have been done, with promising results. 

The observation that the acriflavine-resistant strain requires a much 
higher concentration of acriflavine than the original strain to produce a given 
level of mutagenic activity 1s of interest in relation to a similar comparison 
made by Demeree and Latarjet.* These authors found that strains B and 
Br, which are, respectively, sensitive and resistant to ultra-violet radi- 
ation,’ give the same yield of induced phage-resistant mutants at a given 
dosage of ultra-violet, regardless of the wide difference in survival. The 
apparent contradiction between these observations may be spurious, how- 
ever. It is possible that the resistance of the acriflavine-resistant strain 
may depend upon a permeability change, or upon some other property 
whereby the amount of acriflavine present inside the cell is reduced. The 
survival level may be a measure of the effective concentration within the 
cell, in which case the correspondence between levels of survival and muta- 
genic effect would not be surprising. 

It should be mentioned that in using this procedure it is necessary to 
adjust conditions so as to obtain survivals in the range indicated (10~* to 
10>"), Survivals outside this range rarely give positive results, as there 
seems to be a delicate interplay between the killing rate and the rate of in- 
duced mutation, For any other compound, the optimal survival range 
must be determined empirically 

Summary.—-(1) A method is described whereby exposure of bacteria to the 
action of acriflavine can be sharply arrested, and its residual activity elim- 
inated by the addition of sodium nucleate, 

2) The mutagenic activity of acriflavine, as measured by its ability 
to induce mutations from phage-sensitivity to phage-resistance in strain 
B/rof £. colt, is confirmed by the use of the nucieate-reversal technique. 

(3) An acriflavine-resistant strain of /. coli, derived from B/r, requires 
30 times as much acriflavine as does Br to give a comparable mutagenic 
effect, as well as a comparable level of survival. 

* This study was conducted in part under a grant from the American Cancer Soci 
ety, through the Committee on Growth of the National Research Council 

t The author acknowledges gratefully the assistance of Miss Jessie Flint and Mr 
frank Kennedy in many of these experiments 
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THE NATURE OF BONE AND PHOSPHATE ROCK 
By S. B. Henpricks anp W. L. HILi 


U.S. Burgau Ixpusrry, Sors, AND AGRICULTURAL ENGINRBERING, 
Mp. 


Communicated by M. A. Tuve, September 20, 1950 ‘ 


Carbonate has never been adequately explained as a constituent of bone, 
dentine, enamel, and the closely related phosphorites. Its failure to act 
quite like CaCOs, upon heating for instance, has, since the time of Hoppe- 
Seyler' led to postulated calcium carbonate phosphate compounds as 
essential constituents of bone or enamel. Existence of optically homog- 
enous apatites such as the minerals francolite’ and collophane containing 
carbonate apparently supports this conclusion and has led to the acceptance 
of an apatite structure containing carbonate. 

Preferential solution by acid of carbonate in bone has often been noted*~* 
and has been interpreted by some workers‘ as evidence for the presence of : 
CaCO;y as a separate phase. Crystalline CaCO;, however, has not been : 
observed in bone by the very sensitive test of its high birefringence. A : 
Geiger-counter x-ray spectrometer, moreover, gives no evidence of its 
presence, while one-third as much calcite as would correspond to the 
analysis can be detected when mixed with hydroxylapatite, or naphtha- 
extracted bone. More incontrovertible evidence is afforded by the work 
of Logan and Taylor® which shows that an apatite phase will precipitate, 
from a solution having the inorganic composition of plasma, in the presence 
of carbonate and steadily withdraw CO,~~ from a solution unsaturated with 
respect to CaCQy. ‘ 

Bone, enamel, dentine and collophane, the last being the apparently 
isotropic component of phosphate rocks, give rather diffuse x-ray powder 
diffraction patterns identical, within their definition, with that of fluor- 
apatite. Francolite, the more highly birefringent material of phos- 
phorites, can be obtained as optically clear single crystals containing more 
than 3% CO»? X-ray goniometer photographs of francolite crystals are 
identical with those of fluorapatite, even in very high orders of interference 
maxima. From the structural point of view it seemed necessary for a 
large number of CO; groups to be present in an apatite structure without 
changing the diffraction pattern, even in minor ways, which is very unlikely. 
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The average inorganic compositions of human bone, and enamel, based 
upon a great number of published analyses, can be represented by the 
following formulas’ in which the number of positive tons is taken as 10.0: 


Enamel PC CO) 0.46 | {OH hier (Hol 


Three quite different suggestions for the structural nature of the compound 
present have been advanced in recent years. Trémel and Moller’ and Bale 
and coworkers* consider the material to be hydroxylapatite, 
(OH), with “adsorbed” carbonate and possibly phosphate. Bale* states 
‘Since the diffraction pattern of tooth and bone substance give no indica- 
tion of a carbonate apatite being present, the carbonate ions are also 
probably most abundant in these regions of structural discontinuity that 
may be considered as marking the transition from one crystalline particle 
to another.” Brasseur and Dallemagne in a series of papers published in 
Belgium during the war which were summarized in Nature in 1946" (also 
reference 4) concluded that bone is “tricalcium phosphate hydrate,” 
Cay( POs H,O)s, with CaCO, present as a separate phase and that enamel 
is in part “carbonate hydroxylapatite.”” They state “We cannot admit, 
therefore, that the tertiary with a Ca/P ratio of 1.94 is hydroxylapatite 
with adsorbed PO, ions.’ Hendricks and Hill in 1942? reached essentially 
these same conclusions except that they considered the carbonate of bone 
to be present in the ‘‘tricalcium phosphate hydrate’’ lattice. 

Properties of the several materials pertinent, without further discussion, 
to the question at issue are: 

1. A single phase having the x-ray diffraction pattern of apatite and 
containing carbonate is present even though bone forms in a system un- 
saturated with respect to calcium carbonate. 

2. In human bone this phase is shown by analysis to be a neutral 
substance. It is more basic, approaching hydroxylapatite, in parts of 
enamel, 

3. A large portion of the carbonate in bone can be preferentially re- 
moved by solution in acid without destruction of the apatite lattice. 
Enamel and some ‘‘carbonate apatites’ dissolve umiformly. 

t. Citrate is present in bone as water insoluble form in amounts by 
volume approaching one-half that of carbonate." 

The first and third of these properties indicate rather clearly that car- 
bonate is on the surface of the bone material. Presence of citrate can also 
be explained in this way instead of being merely “adventitious.”” Bone 
also contains magnesium and sodium that might be expected to replace 
Ca** in the apatite structure on the basis of ionic radii for Na* and 
Mg**. These two are preferentially dissolved by acid,® however, indicat- 
ing their location on the lattice surface. The difficulties of charge with 
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Na* replacing Ca** and essential absence in solutions near neutrality of 
magnesium as Mg**, rather being Mg(OH)*, accordingly do not arise. 
There is still a size limitation for the surface, however, since Na* is taken 
up preferential to K*. 

The nature of ‘carbonate apatites’ has been most confusing in its 
bearing on the type of compound present in bone. Lack of preferential 
solution in acid of carbonate in enamel was essentially the reason for 
Brasseur and Dallemagne considering the compound to be ‘“‘carbonate 
hydroxylapatite.” Another possibihty, however, is that the carbonate 
groups of enamel and francolite and the citrate of enamel are also present 


4 
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FIGURE 1 


Entrapped surfaces in ‘carbonate apatites” as present in enamel 
and francolite, schematic. 


on surfaces which are entrapped as crystallites grow and are surrounded by 
phosphate such that the final material has many occluded surfaces as 
schematically shown in figure 1. In the case of francolite, and to a lesser 
extent in enamel, the apatite lattice can be continuous throughout large 
crystals. The high density of francolite’ indicates that voids cannot 
constitute more than 1°), of the total volume. 

Evidence in regard to entrapped surfaces for francolite and collophane 
comes from some peculiarities of their fluorine contents which have not 
previously been explained. These substances® '* contain fluorine in excess 
of that required for an apatite structure, sometimes by a factor as great as 
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1.0, while fluorapatite has not been found to contain excess fluorine. 
Now if the surface occurs at primitive translations and all fluorine positions 


FIGURE 2 


Excess fluorine in an apatite-type crystal due to lattice limitation 


are filled then part of the surface flu 


orime atoms would be in excess of the 


unit cell requirements as shown in figure 2. Francolite is always observed 
to be secondary in formation to collophane and as it forms with a decrease 


TABLE 1 


Errecr or Fracmenr Size on Surrace ArRBA AND SoLuTION rn Acip or Some Car- 
BONATE CONTAINING APpaTITE SUBSTANCES 


PARTIAL SOLUBILITY IN ACID, 


FINENESS SURFACE PER OF TOTALS-———. 
Steamed bene 100-200 66.2 18.6 6.5 6.5 
39 66.2 15.5 5.9 5.4 
63.7 

Francolite Relow 6 micron 7.3 26.1 7.2 2.9 
100-200 3.2 10.5 6.6 3.9 
20- 35 1.6 6.5 3:7 3.5 

Knamel Through 325 6.8 

200-325 3.7 

120-200 3.2 

60-120 2.1 


* Analyses by J. H. Caro. 
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Other evidence for entrapped surfaces can be obtained from surface area 
measurements by nitrogen absorption (B.E.T.). These areas should de- 
pend markedly upon screen mesh size of material under test if entrapped 
surfaces are uncovered by breakage as indicated in figure 1. That this is 
true for enamel] and francolite, but markedly less so for bone as required by 
the preferential solution of its carbonate in acid is shown by results in table 
1. These surface areas were measured in the laboratory of Dr. V. R. 
Dietz of the National Bureau of Standards who, independently, had made 
the observations on enamel for Dr. R. C. Likens of the National Institutes 
of Health. 

Attention is now turned to some detailed consideration of the surfaces. 
Bone, after removal of most of its organic matter, has a maximum surface 
of about 80-100 m.*/g. (Reference 15 and extensive unpublished informa- 
tion of V. R. Dietz for bone ash) and contains about 7.5 X 10” CO, 
groups/g.’. The corresponding area per group, assuming a surface area 
of 100 m.*/g., the larger value being taken in allowance for entrapped 
surfaces, is 13.3 A* which is about the value to be expected if the carbonate 
(and possibly citrate ions) forms a single surface layer containing cations 
and the maximum surface is realized. 

Now consider the francolite sample analyzed by Gruner and McConnell® 
which contains 3.360% CO, and 1.08 times more fluorine than required for 
apatite. The excess surface fluorine is given by 2x 4+- 1 = 0,08x*, where 
x is the number of primitive translations in a ‘‘crystallite’’ which is assumed 
to be equal in the x and y directions, figure 2. Thus x is about 25 corre- 
sponding to 234 A (25 X 9.35 A). The calculated ‘surface area’ /g., made 
up chiefly of that entrapped, is about 62 m.*/g. The “external surface”’ 
is surely less than 1 m.*/g. A completely independent estimate of surface 
area is afforded by the carbonate content. For 3.36% CO, and an area of 
13.3 A* per CO; group this is 61 m.?/g. While the agreement with the other 
value is perfect the combination of the methods can only be trusted to 
about 25°), the greatest uncertainty being about the surface area per 
COs group. 

The nature of “tricalcium phosphate hydrate’’ having an apatite diffrac- 
tion pattern, and serving as a prototype for bone, is now considered. This 
material can be formed by precipitating calcium phosphate in a system 
having a pH near that of plasma or by hydrolysis of CaHPO,, i.e., under 
conditions where H PO,~~ is a predominant phosphate ion in solution. The 
composition of this material determined by independent groups is Cag- 
or Cag( )e(H20);.5," the question of the exact water 
content of course being difficult to settle. A preparation made by slow 
hydrolysis of CaH PO, in a system more acid than pH 5.0 contained crystals 
measuring about 2 microns in length which were shown by electron micro- 
graphs to be well formed lath-like prisms with pyramidal faces. A ma- 
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terial of this type gradually lost 4.70°) HyO (calculated for Cag{ POy)s( HzO): 
3.73°) HyO) between 100° and 800°. Above SOO° the substance is essen- 
tially anhydrous and gives the diffraction pattern of the beta modification 
of Cag( PO, (stable below 1115°C. and possibly a stable phase in solution) 
which verifies the molal ratio. Hydroxylapatite behaves in a markedly 
different way in retaining its full contént of 1.799% H,O upon heating at 
1000°. It is dehydrated fully only above 1300°. 

Is this substance really what it appears to be or is it hydroxylapatite 
as was suggested many years ago by Trémel* (also references 6 and 9) 


with surface material modifying the composition? The “formula” on 
this basis would be 


Carol IH)» Surface ( P( 
Surface (H )e.00 


and would have to be approximated, on the basis of present evidence, 
over a pH range from about 4.6 to 7.3 in the precipitating system. The 
particular composition might vary somewhat and would be determined 
possibly by the surface and the pH. Surface areas of several materials 
made by precipitation varied from 26 to 64 m.*/g. With surfaces, sepa- 
rated by perhaps ten primitive translations and entrapped in part, covered 
with HPO,~~ ions, behavior upon heating might parallel that observed; 
namely, gradual loss of water in a manner very different from that of 
hydroxylapatite. In the case of bone, HPO,~~ would have to be present 
on the surface together with carbonate and citrate. Experimental results 
at hand cannot safely resolve the question, although the authors, as well 
as Brasseur and Dallemagne, have previously considered the evidence 
adequate for “tricalcium phosphate hydrate.’’ The observation that the 
poorly developed x-ray powder diffraction pattern is closely the same as 
that of hydroxyl apatite should surely be given full weight,*® but it alone is 
not conclusive. 

This surface chemistry of bone and related materials markedly affects 
composition and is a result or probable cause of the large surface area. 
It can be generalized. An insoluble calcium compound can form surface 
compounds with other ions having insoluble or poorly ionized calcium salts 
(e.g., citrate), Thus CaCO, would equally be expected to form surface 
phosphates and citrates in much the same way as calcium phosphates 
form surface carbonates. Such reactions which are important for phos- 
phate fertility of caleareous soils have been observed in tracer experiments 
by Dr. Sterling Olson of this Bureau. 

The phenomenon discussed here, which can be studied in many ways, 
has broad implications in the formation and properties of bone and teeth. 
\n explanation is afforded for the unusual relationship of fluorite, collo- 
phane and francolite in phosphate rock. Interpretation can be given to 
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occurrence of uranium and other groups as surface material in phos- 
phorites."" An understanding is afforded for the preferential use of certain 
phosphate rocks, particularly those from North Africa having high car- 
bonate contents, for direct application to some crops." 

Modification of composition by small particle size has also been found for 
hydroxylapatite by P. W. Arnold of Cambridge University (personal 
communication). He observed that precipitates of hydroxylapatite con- 
tain (OH)>~ in excess of the formula requirements, such as to yield more 
than 2.0°% water between 900° and 1400°C. (calculated 1.79% HyO). 

Summary.—The inorganic compounds of bone have a large surface on 
which carbonate, citrate, magnesium and sodium are located. Carbonate 
is present on entrapped surfaces and not as constituent within the lattice 
of the mineral francolite, the typical “carbonate apatite,” that also by 
lattice limitation, contains fluorine in excess of the fluorapatite require- 
ment. Enamel is similar to francolite in having entrapped surfaces. The 
over-all composition, including lattice and surface, of bone is that of a 
neutral calcium phosphate while enamel is more basic, approximating 
hydroxylapatite. The compound forming the lattice of bone remains to 
be established but it might well be hydroxylapatite with an excess of 
phosphate on the surface. 
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LIMITATION OF ELECTRON DENSITY BY THE PATTERSON 
FUNCTION 
By M. J. Buercer 
Lasoratory, Massacnuserrs insrirurs or TecHNnoLocy 
Communicated by W. J. Mead, November 1, 1950 


In several recent papers'~* it has been shown that the vector set can be 
solved for its fundamental set, and this implies, in a broad way, that a 
crystal structure can be solved from the intensities of its x-ray diffraction 
spectra. In particular, it implies that hitherto unsuspected relations 
exist between the Patterson function and the electron density function. 

One such relation has been demonstrated: If an image of any polygon, 
whose vertices (n 2) have coordinates xyy2,, . . ., Can be 
found, one can construct the function 


= P(xyys:) K ..., (1) 


where P(x,y,%), etc., are Patterson functions. Then if there are no equal 
vectors between points in the fundamental set, there are no fortuitously 
coincident images in the vector set, and the following relation holds,’ 
except for an exaggeration due to the selected image point: 


p(xy2) = K VII, (xvs), (2) 


where A is a constant defined by the chosen image. 

For actual crystals, the electron density varies continuously with the 
coordinates, so that the Patterson synthesis contains a continuous over- 
lapping of images. These arise not only from Patterson vectors radiating 
from a particular origin, but also from Patterson vectors reaching into a 
particular cell from the origins of the neighboring cells. The Patterson 
function for an actual crystal therefore contains extensively overlapped 
and coincident images. This feature always degrades the information 
obtainable from (1) in a specific way, namely by increasing the values of 
one or more P’s and therefore increasing the value of I]. For an actual 
crystal, therefore, an inequality corresponding to (2) holds, namely 


p(xyz) KVII,(xys). (3) 


In spite of this unfortunate feature, the IT function is useful in the 
interpretation of Patterson maps. Only high ground can be occupied by 
atoms, and low ground can be regarded as unpermitted territory. In 
instances where the overlapping is minimized, the function is most useful. 
For example, a collection of heavy atoms embedded in a matrix of lighter 
ones can be readily located by the II function. The reason for its power 
in such an instance is that the characteristics of the set of heavy atoms 
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approaches that of a set of pure points reasonably separated in empty 
space. 

A particularly bad feature of the II function depends on the fact that 
a Patterson map usually has a more or less continuous non-zero background 
which arises from at least partial overlapping of the images of atoms. 
When a line image, for example, encounters a multiple peak at one end, 
it usually has this background at its other end, so that the product has some 
moderate value. Because of this, the II function often produces a ghost 
when one point of the polygon passes over a large Patterson peak, in spite 
of the fact that the polygon location is not an image position. Due to the 
presence of such ghosts, the II function usually has more peaks than the 
electron density function. 


FIGURE |! FIGURE 2 
Centrosymimet- Vector set of the fundamental 
rical fundamental set of Fig. 1, showing weighting 
set. of points (circles) and also 
values of the M; function 
(crosses). 

Obviously, the II function is not the only image-seeking function which 
can be devised to locate atoms. Any function which drops to zero when 
the weighting on at least one of the polygon vertices is zero can be expected 
to map the electron density with more or less faithfulness. A particular 
function which avoids the disadvantage of the II function mentioned in the 
last paragraph may be called an VW function. To form the M, function, 
one records at xyz the minimum value of the » values of the Patterson 
function which occur at the n vertices of the roving image polygon. -It 
will be demonstrated that this function has the interesting property that 
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pixys) AM,(xys), (4) 


using the same notation as for (2). Inequality (4) is more powerful than 
(3) 

The general properties of the M function can be readily appreciated if 
it is applied to the case of a polygon image of nm points which are related 
by symmetry in the crystal. For simplicity and clearness, consider n = 2, 
that is, a line image, and apply this image to a vector point-set. Let the 
point at the end of the origin line-image be a product (“‘interaction’’) 
between centrosymmetrically related points in the crystal. This is illus- 
trated in figures land 2. Figure | shows the points in a cell of the crystal 
and their weights; figure 2 shows the vector points arising from the points 
within the single cell shown in figure 1. The dotted lines of figure 2 show 
the line images, and the cross on each one shows the location of a peak of 
the AM, function. The labeling indicates the weights of the peaks.  Inci- 
dentally, note that all points of the vector set which are “interactions” 
between centrosymmetrical points of the crystal set have unit weight, 
while other “interaction’’ points have double weight. It will also be 
observed that the locations of the peaks of the M/, function are the same 
as the points of the crystal set, and their magnitudes are the same except 
for a factor 2a, which depends on the weight (corresponding to electron 
density) of point a, which was chosen for the first image point. An excep- 
tion in weighting in the 47 function occurs for the point reproduced by the 
line image containing the origin; this has unit weight. The exception 
occurs only in the case of the point associated with the vector peak chosen 
for the line image to be used for the decomposition of the vector map. 

The M function not only has the advantage mentioned above over the 
II function, but also a related advantage of usually giving a more correct 
value when the image polygon comes to an image location for which one of 
the vertices lies on a multiple peak. As the Tl function ranges over such an 
overlapped peak, it lays down an exaggerated peak, but the 7 function 
lays down a peak of normal magnitude provided that at least one vertex 
of the image polygon is not overlapped. Under such conditions, therefore 
(except for the neighborhood of the atom a used as the initial image 
points) the W function has the important property 


p(xys) = ~~ M(xys). (5) 


In the case of an actual crystal, the fundamental set becomes the electron 
density map, and the vector set becomes the Patterson map.’ There is 
usually overlapping to some extent at all vertices of the image polygon in 
the Patterson map. But even when this occurs, it merely increases the 
value of M, and the inequality 
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l 
pixys) < (6) 


=Pe 
is valid. Relation (6) implies that the Patterson function places an upper 
allowable limit of the electron density at every point in the cell. Since 
the electron density at the center of an atom is a characteristic of an atom, 
this relation will permit locating atoms in instances where the Patterson 
is not complicated by too much overlapping. In any case, it permits 
eliminating regions where a specific atom cannot be located, on the ground 
that the electron density at the atom center is greater than that of the MW 
function in the region. 

In crystals having symmetry other than an inversion center, a centro- 
symmetrical fraction of the cell is repeated by other symmetry elements to 
form the entire cell. When the image polygon is chosen so that it does not 
conform to the Patterson symmetry at the origin, it does not, in general, 
encounter symmetrical situations at equivalent locations in these fractions 
of the cell. Nevertheless, if the Patterson map contains no fortuitous 
overlaps, any image function synthesizes the true symmetry of the crystal. 
But if there are fortuitous overlaps, the image function may synthesize 
a degradation of the true symmetry. For example, consider a two- 
dimensional Patterson map of symmetry (2mm, and let a line image be 
chosen which, at the origin, makes an oblique angle with both symmetry 
lines. In reproducing symmetrically located peaks on the two sides of a 
symmetry line, the sloping image line uses two peaks on one side of the 
symmetry line and two peaks on the other side of the line, respectively. 
These two pairs of peaks are not symmetrically related unless the line 
image conforms to the symmetry of the Patterson map. Since, in general, 
each pair of peaks is involved in different overlaps, any line function re- 
produces the symmetry of the crystal as degraded by the overlap charac- 
teristic of the cell fraction. In such instances, the minimum value of 
either the Il or the M function at equivalent points in different cell fractions 
can be combined to give a better minimum value of Il or WV, thus increasing 
the power of the particular inequality (3) or (6). 

Since (6) is quantitative, it requires the Patterson synthesis upon which 
it is based to be quantitative. Any feature which reduces the accuracy 
of the Patterson synthesis therefore reduces the effectiveness of (6). A 
particularly bad effect has been noticed which is caused by the limited 
number of experimental #’'s available and the consequent termination 
of the Fourier series used for the synthesis of the Patterson map. This 
gives rise to false sets of undulations corresponding to optical diffraction 
fringes. A set of such undulations surrounds each strong Patterson peak, 
but the undulations are most noticeable surrounding the great origin 
peak. The ring-shaped trough and the following ring-shaped hill which 
immediately surround the origin are very objectionable and tend to distort 
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the value of the VW function as the image polygon passes over them. It is 
recommended that before computing an M function from Patterson data 
these undulations be eliminated along the lines suggested by James.‘ 

As described above, the A/ function was applied to polygon images 
composed of points due to symmetrical ‘interactions.’ It can also be 
applied to image polygons containing points due to non-symmetrical 
“interactions.” The selected image is then a polygon whose n points 
have weights, say, aa, ad, ac, . Toform the M function for this polygon, 
one first prepares m Patterson maps of the sampling of the Patterson 
function sealed, respectively, by factors 1/a, 1/6, l/c, ..., and lets one 
point of the imaye polygon rove over the corresponding map. When the 
polygon is located at xys, the minimum value which occurs at the n, points 
is then recorded at xyz, and the following relation holds 


p(xys) Mixys). (7) 


Buerger, M. J., Sets," Ada Cryst., 3, 87-97 (1950). 

* Buerger, M. J., “The Crystallographic Symmetries Determinable by X-ray Diffrac- 
tion,” Proc. Nati. Acap. 36, 324-329 (1950). 

* Buerger, M. J., “Some New Functions of Interest in X-ray Crystallography,” 
Ibid., 36, 376-382 (1950) 

* James, R. W., ‘False Detail in Three-Dimensional Fourier Representations of 
Crystal Structures,” Acta Cryst, 1, 132-134 (1948) 


SOME LINEAR OPERATORS IN THE THEORY OF PARTIAL 
DIFFERENTIAL EQUATIONS* 
By S. BERGMAN AND M. SCHIFFER 
HARVARD UNIVERSITY 
Communicated by Marston Morse, September 29, 1950 
1. Let D, be a finite domain in the x, y-plane bounded by an analytic 
curve Cy. Let D, be a proper subdomain of D, with analytic boundary 
curve C;. Denote finally the difference domain bounded by C, and C; 
by Dy 
Consider the partial differential equation of elliptic type 


Au = g(P)u, (1) 


where g(/) is a positive continuously differentiable function of the point 
P we (x, y) in Dy + Cy. We introduce the scalar products 


Ejju,v} = SS (grad u-grad + queojdr, = 0,1,2, (2) 
D, 
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dr = area element in D, 


and the three Hilbert spaces =, consisting of all solutions « of (1) with a 
finite norm E(u) = E,{u, wu}. Let finally N,(?, Q) and G,(P, Q) denote 
Neumann's and Green's functions with respect to the differential equation 
(1) and the domain D,. 

We define the kernels 


Q) = N&P, Q) — GAP, Q), LAP, Q) = NP, Q) + Q) 
(3) 


which play an important role in the theory of the spaces Z, [1,a@ — d; 2] 
and each of which is defined in the domain D,. We may construct now 
discontinuous solutions of (1) which are defined in D,, composed in a simple 
way of the kernels (3) and possess various interesting properties. We set 


L,(P, Q) — LA(P, Q) for P, Oe D, 
= — LP, Q) for Pe D;, Qe Do or P€ Do, Qe D, 
LAP, Q) L.(P, for P, € Ds. 
(4) 
Similarly, we define 
K,(P, Q) KAP, Q) for Qe D, 
k(P, Q) = Ki(P,Q) forPeD,, QeDoor Pe Do, QeD, 
Kol P,Q) — Kil P, Q) for P, Qe De 
(5) 
Using the characteristic properties of the various Green’s and Neu- 
mann’s functions, the following identity can be established: 


E,{l(P, Q), R)} = R(Q, R), Q,ReD,but notonC, (6) 


Here the convention will be used that if kernels with various argument 
points are multiplied the operation shall always be understood with re- 
spect to the common argument letter. Similarly, we can show 


E,{l(P, Q), R(P, R)} = R), Q, R « D; but not on C,, (6’) 


E,{k(P, Q), R(P, R)} = k(Q, R), Q, Re Dz but not on C, (6”) 


2. Let 2 denote the Hilbert space of all functions in D, which represent 
in Dy an element of 2» and in D, an element of Z;. The metric in = will be 
based on the scalar product F,{u, v}. 

We define in © the linear operators [cf. references 1d and 3} 


T,(P) Q), u(Q)}, S,(P) E,{K(P, Q), u(Q)}, (7) 
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which transform the space = into itself. They satisfy the following rela- 
tions 


Ex{u,v) = ExlT., Te} + Esl So} (8) 
S| = 0. (8’) 


These results suggest the imtroduction of the two linear operators in =: 
= T, + S, = T, 3, (9) 

which have the property of preserving the metric 
S = vf, S| = (10) 


The operators S,* and @,~ have further the following involutory char- 
acter 

Ife = Sy, thenu = S,* andifw = S,~ thenu = S,~. This result 
shows that ©* and @~ are one-one mappings of E onto itself. 

ach funetion «a ¢ 2 can be decomposed in a unique way into 


uiP) = a(P) 4+ b(P) (11) 


where « and a(P) is orthogonal to The operators T and S per- 
form this orthogonal decomposition of Z. 

The operators S*, ‘S~ can easily be interpreted as follows: Given an ele- 
ment « « Z, which is continuously differentiable in the closures of Dp and 
D,, S,* will be that function of 2 which has on C, the same boundary values 
as u, which has along C; the same discontinuity a su but whose normal deriv- 

ative has the discontinuity of opposite sign on C, than has ~ : 

Similarly S,~ has on C, the same normal derivative as u, its normal 
derivative has the same jump through C, as that of «, but its discontinuity 
across C, is at each point minus the jump of « over that line. 


Let we be identically zero in Do. In this case = and 
we derive from (8): 


2 Ey(T,). (12) 

This inequality holds for every function « «,. We derive from it, in par- 
ticular, that the eigenvalues of the integro-differential equation 

uP) = Q), Ve D, (13) 
are ull greater or equal to | in absolute value. A finer study of the condi- 
tions for the equality sign shows that always 
> 1 (14) 
must hold. This result plays an important role in the theory of the kernels 
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4. Let us define the linear operator in 2, 
s,*(P) = Q) + Kal P, Q), u(Q)}. (15) 
From the continuity properties of the kernels / and A; one can easily show 
that this operator is completely continuous and, hence, treat the eigenvalue 
problem 
= «, Sy *(P), PeD, (16) 


by the classical methods of Hilbert space theory. We obtain a discrete 
spectrum of eigenvalues «, and corresponding eigenfunctions u,(P) which 
we may assume to form a complete orthonormal set in =;. By the reason- 
ing of the last section, one can easily prove 


> 1. (16’) 
We have in D, the Fourier developments 
a l 
Ki(P, Q) = u,(P)u,(Q), KP, Q) + Q) = u,(P)u,(Q). 
p= l 
(17) 


Using the definitions (3), we easily derive from (17): 


l 
GAP, Q) — GAP, Q) = (: ) u,(P)u,(Q). (18) 


pel 


Consider the @*-transforms of the eigenfunctions u,(P) ¢ =, which are 
extended into functions of the class = by defining them to be identically 
zero in Do. Let 


v,(P) = S.,*(P) = —2E,1G,(P, Q),u,(Q)}, (19) 

The {v,(P)} form an orthogonal set of functions in Xo; they satisfy 
v,(P) = Q) + K:(P, Q), »,(Q)} (20) 
and thus play a very similar role in Xo as did the u,in Z,. We calculate 


Eolv,, ve) = (1 (21) 


All v,(P) vanish on C,; they form a complete set of functions in the sub- 
space of X» which consists of all elements vanishing on C,. The function 
GP, Q) — Gol P, Q) belongs to this subspace and can, therefore, be de- 
veloped into a Fourier series in terms of the v,(P): 


2 Ky 


Since the 9,(P) can be calculated from the u,(P) by elementary integra- 
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tion if Gy is known, we see that the u,(P) permit a representation of the 
Green's functions of the complementary domains D, and D,. 
Using instead of s* the linear operator 


s, (P) = Q) Ki(P, Q), u(Q)} (23) 
we can develop an analogous theory for Neumann's functions. 


5. The function w,(?) «2 which has the value ~ u,(P) in D, and v,(P) 
Ky 


in Dy is the @*-transform of the function u, « 2 which has been extended 
from 2, by defining it identically zero in Dy. Using the known discontinui- 
ties of the function u« across C, we derive 


v,(P) = (1 ) ar), PeC, (24) 
K, 
ov,( P) ( ‘) ou,(P) 
= {1 PeC, (24’ 


Thus, the boundary values and the values of the normal derivatives of the 
functions v, ¢ X» and uw, « X, are proportional along the common boundary 
C, of their domains of definition. 


The function w,(P) ¢ > which is defined as (; _ i) u,(P) in D, and 


as 0,(P) in Dy has the following property: It is continuous in D,, vanishes 
on the boundary C, of D, and satisfies along the curve C, the discontinuity 
condition 


ov" x, | 


Thus, the eigenvalue «, can also be characterized by an eigenvalue prob- 
lem for discontinuous solutions of (1) of the Stekloff type [4]. 


* Work done under Navy Contract N5ori 76/XVI NR 043 046, at Harvard Univer 
sity 

' Bergman, S., and Schiffer, M., (a) “A Representation of Green's and Neumann's 
Functions in the Theory of Partial Differential Equations of Second Order,’ Duke 
Vath. 14, 808-038 (1947). “On Green's amd Neumann’s Functions in the Theory 
of Partial Differential Equations,” Bull, Am. Math. Soc, $3, 1141-1151 (1947). 
“Kernel Punctions in the Theory of Partiuil Differential Equations of Elliptic Type,” 
Duke Math J., 1S, 585-566 (1948) (d) “Kernel Functions and Conformal Mapping,” 
Compositio Math., 8, (1950) 

* Schiffer, M., “The Kernel Function of an Orthonormal System,” Duke Math. J., 13, 
A240 (1046 

* Schiffer, M., and Spencer, D. C., 'Functionals of Finite Riemann Surfaces,” Prince- 


ton University Press, to appear 
Stekioff, W., "Sur la Théorie des Fenetions Fondamentales,”’ C. R. Acad. Se. Paris, 
128, 
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GENERALIZED GROUP ALGEBRAS 
By R. CrercuTron Buck 
DepaRTMENT OF Maruematics, Tue University oF WISCONSIN 


Communicated by Saunders MacLane, November 2, 1950 


The purpose of this note is to describe a natural banach algebra asso- 
ciated with any topological group; the theory of this algebra and its 
subalgebras subsumes much of the standard L' theory, and is independent 
of the existence of a measure. Proofs will appear elsewhere.' 

Let G be an arbitrary topological group, not necessarily locally compact, 
and let K be the complex field. Let C be the banach space of all bounded 
uniformly continuous functions on G to K, with norm || sup] fix)|. 


Let B(C) be the banach algebra of all bownded linear operators on C. For 
b eG, let U, be the left translation operator which sends a function f into 
» whose value at x is f(bx), and let U® be the corresponding right transla- 
tion. By our choice of C, U, and U* are in B(C) and have norm 1; the 
mapping x — LU, (x — U*) is an anti-isomorphism (isomorphism) of G 
in B(C), not in general continuous. 

Derinirion: Let B be the set of all T in B(C) which commute with ail the 
operators 

B is a closed subalgebra of B(C); it is the centralizer in B(C) of the 
subalgebra generated by the left representation of G. Since U* and U, 
commute, 8 contains the nght representation of G in B(C), and we may 
therefore regard G as embedded in 8. Consider G, the set of all bounded 
continuous homomorphisms of G into A°, the multiplicative group of K. 
Clearly, G is a subset of the unit sphere in C. 

THEoreM 1. A function f in C is a character or a multiple of a character 
if and only if f is simultaneously an eigenfunction for every operator in B. 

Thus, B determines G and hence G, when G is such that the duality 
theorem holds. Regarding G as embedded in ¥, it has two topologies, 
the group topology carried over by the embedding, and the relative norm 
topology which it acquires as a subset of B; the latter is usually discrete. 

TuroreM 2. Any character of G has an extension to a ring homomorphism 
of B into K, continuous in the norm topology. 

li heG, then the extension H is defined by H(T) = A, where A is the 
eigenvalue of h under T. The kernel of H is simply the set of all T in 
B with 7(k) = 0, the origin of C. 

Corottary: G has a natural one-to-one embedding in Hom(®:K) the 
set of all continuous ring homomorphisms of B onto K. 

Clearly, any homomorphism H of 8 onto K has a contraction to G which 
is a group homomorphism of G into K°. (This remark, together with the 
preceding theorem, essentially contains the results obtained by Kakutani 
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and Kodaira.") However, this homomorphism is not necessarily in G 
since it need not be continucus in the group topology of G. Under certain 
conditions, this continuity can be forced. For example, this is the case 
if there is an operator 7) ¢ B such that /(7,) = 1 while the map x — UT, 
is continuous. This is the mechanism underlying the usual! L' theorv. 

Since B is large, the study of its subalgebras is of importance. Let 
YW, be the closed subalgebra generated by the image of G in its natural 
embedding. We call a subalgebra WM total if T(f) = 0 for every 7 ¢% 
implies that f = 0. It is convenient to single out the class of linear func- 
tionals L on B* having the form L(T) = Tig)(e) where ¢ « C and ¢ is the 
neutral element of G.* We term such regular functionals; they form a 
subspace of the dual of B*, usually not closed. The homomorphisms which 
anse as extensions of those in G are regular. Given a subalgebra MA ¢ B, 
G has a natural one-to-one embedding in Hom(%:X) which in certain 
cases 18 onto 

Turorem 3. Let be a closed total subalgebra of such that (21) AM A 
and WNy ¢ A, (2) the map x — U'T ts continuous for each T eA, (3) if H 
is any homomorphism of A onto K, then H(U*°T) = H(TU") for all beG 
and T ¢ M, (4) every member of Hom(A:K) is the limit of regular functionals. 
Then, every homomorphism of % onto K arises as the unique extension of a 
member of G 

lo show the connection of the preceding with the usual theory, let ¥ 
be the dual space of C. The correspondence T — F, where F( f) = 7T(f)(e) 
may be shown to be a norm preserving isomorphism of the banach spaces 
Wand ¥ Since B has the additional structure of an algebra, % can be 
given the same structure, Denoting the resulting multiplication for 
functionals by *, F,* #; turns out to be the customary convolution of the 
functionals ‘and Fy, in the sense of Radoti measures.‘ The operator 
L” in B corresponds to the point functional 6, where 6,(f) = f(b). The 
theory of the Fourier transform assumes the following simple form: given 
Fe, let P be its restriction to the subset G ¢ C. The importance of 
discussing the span of G in Cis evident. For the usual L' theory, choose G 
to be locally compact abelian, and let YM be the set of all operators of the form 
gix) = Sflxy)e(y) dy where ge L(G) and feC. The hypotheses of 
Theorem 3 are easily verified. 

The norm topology on C is not the only usable one; in particular, that 
of uniform convergence on compact sets deserves special consideration. 
The results for this case are similar, and will appear elsewhere. 

' There seems to be some relation between the present approach and that in a recent 
abstract of Hewitt and Zuckerman, Auli. Am Math. Sec., $5, 1057 (1949). See also 
1. B. Segal, Ann. Math, $1, 203-208 (1950). The two-sided regular representation of 
a unimodular locally compact group 


* Kakutani and Kodaira. Normed ring of a locally compact abelian group, Proc. 
Imp. Acad. (Tokyo), 19, 360-365 (1943). 
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* B* is the underlying banach space of the algebra B. 
* Weil, L’integration dans les groupes topologiques, Actualitiés Sci. et Indust., Paris, 
1940. 


LIMIT CYCLES OF SYSTEMS OF THE SECOND ORDER 
By G. F. D, Durr 
Princeton Untversiry, Princeton, N, J. 


Communicated by S. Lefschetz, October 2, 1950 


The present paper deals with a number of questions about the limit cycles 

of systems 
& = P(x, y), ¥ = Qt, »). (1) 

The importance in non-linear mechanics of these limit cycles, or isolated 
periodic solutions, is well known. It has led to a variety of methods for 
finding these solutions which yield some information in noteworthy special 
cases. In Princeton lectures in the fall of 1949, Lefschetz introduced and 
applied to certain questions the following method: a rotation by a fixed 
angle a being applied to the vector field (Y, Y) in the phase plane of the 
system (1) gave rise to a new vector field (Y cos a — FY sin a, X sin a + 
Y cos a), whose trajectories were then compared to those of the original 


system. A method essentially like that of Lefschetz is utilized in this note, 

and has made it possible to answer systematically many questions con- 

cerning limit cycles of second order systems. ; 
Consider then a family of systems 


% = P(x, y, a), y = Q(x, y, a). (2) 


We shall say that these equations form a complete family, and that the 
associated vector fields in the (x, y) plane form a complete family of rotated 
vector fields in the plane, if the following conditions are satisfied : 

(a) Pand Q satisfy Lipschitz conditions with respect to x and y, and are 
differentiable functions of a. 

(b) The zeros of P* + Q* remain fixed in the (x, y) plane as a varies. 

(c) The condition 


is satisfied at all points where P, + Q* > 0. 
(d) P and Q are periodic functions of a with minimum period 2”, and 
furthermore 
P(x, y,a + 4) = —P(x, y, a), 
Q(x, y,a + = —Q(x, y, a). 
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As @ varies, one may picture the field vectors as turning, at a positive rate, 
all in the same sense, and in such a way that, in the interval 0 < a < #, 
all have turned through the angle x. Thus the original field, but with the 
vectors reversed, appears again. An example is the set of uniform rota- 
tions of the field vectors by an angle a. The critical points remain fixed, 
according to condition (6). We shall study the phase portraits of these 
fields, and will show that they are closely related. 

We note that the index (in the sense of Poincaré) of any closed curve in 
the plane, with respect to any field of a complete family, is independent of 
a. It may also be verified that if a change of the dependent variables of 
the form x, = x(x, ¥), ¥ = ¥i(x, y), with functional determinant every- 
where different from zero, is made, then the differential equations in x, and 
y, also constitute a complete family. 

The following theorems give fundamental properties of complete families 
as defined above: 

Turorem |. Limit cycles corresponding to different values of a do not 
intersect. 

Tueorem Il. Lamat cycles which are completely stable or completely un- 
stable expand or contract monotonely as a varies ina fixed sense. Semistable 
cycles either split into two diverging cycles or disappear. 

Tureorem III, Hach elementary critical point of positive determinant A 
either generates or absorbs exactly one cycle in the rangeO [ a < #. 

We sketch a proof of Theorem IIT for the case in which the rotations are 
uniform over the plane. Let 


x= X cosa — ¥ sin a, 
y= sina + cos a, 
be the equations of the family, where 
X = ax + by + X,, 
= ox + dy + 


and X; and }, contain only terms of the second and higher orders in x and 
y. The characteristic roots of the linear equations of first approximation 
near the origin are given by 


Mt — [(a + d) cos a + (b — c) sin a] A + ad — be = 


Writing 
A = ad — be, 


= (a + d)*? + — c)’, 


a +d 


tané = — ‘ 
c—b 
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we can write the characteristic roots in the form 


Vi int 8) 


9 


A= + > sin (a — 3) = 


If A > 0, corresponding to an index of +-1, these roots assume both real 
and complex values, since4A <I*. For a = 6, the roots are purely imagi- 
nary, and the originisacenter. For nearby values of a, the origin is a focus, 
or spiral point, with the stability of the spirals changing sense when a = 6. 
According to Theorem I, the generation of such a cycle can take place only 
once in a half revolution of the field. 

With the aid of these theorems, it can be shown that regions of certain 
types in the phase plane are swept out by the cycles of a given complete 


’ family. Examples of such regions are: the interior of a given cycle con- 


taining a single critical point, which is elementary, or an annulus free of 
critical points which is bounded by two cycles of the family not necessarily 
belonging to the same value of a. In some cases the whole plane is covered. 
These results, together with Theorem I, lead to criteria for the non-existence 
in such regions of limit cycles of other equations of the given complete 
family. 

The theorems quoted above also suggest the following method of locating 
limit cycles. Embed the given equation in a complete family of rotations 
and calculate the values of 4 for the elementary critical points of index 
+1. The expansion of a cycle generated at one of these points may be cal- 
culated as follows: When the cycle is very small, the linear terms in the 
expansions of Theorem III predominate in the right-hand sides of the 
equations and it is possible to use Poincaré’s method of generating cycles 
to determine the size. This method picks out the closed solution of an 
equation with small non-linear terms from among the closed solutions of 
the corresponding linear equation. The rate of expansion may subse- 


d 
quently be found with the following formula for . the rate of expansion of 
a 


the cycle along its normal: Let T be the period and A the exponent of the 
cycle. Let v(t) denote the “velocity” (P? + Q*) ve Then, 


t+r og oP | 
da v(t) [1 — exp(—AT)] 


The cycle may disappear in one of the following ways: (a) by coalescing 
with an exterior or interior cycle of opposite stability (in this case, reverse 
the variation of a and follow the other cycle); (+) by tending to a saddle 
point or other critical point; (c) by tending to infinity in one or more di- 
rections. 
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If the cycle in question fails to persist until the value of a correspond- 
ing to the given equations is reached, then Theorem | asserts that no part 
of any cycle of this system lies in the region swept out during the expansion. 
For any value of a the cycle may become a band of positive width of closed 
cycles. When a has run through an interval of length +, each cycle must 
either have disappeared or else have taken the place of another cycle of 
the orginal field. Thus the method will always find those limit cycles 
which enclose one critical point with A > 0, and it provides a criterion of 
non-existence, in case there are no such cycles. 

In more complicated situations where the cycle might surround several 
critical points one must establish the existence of a cycle of the family for 
some value of a. This may be accomplished by finding a closed curve with- 
out contact with the given field (a Bendixson curve), and then choosing the 
complete family in such a way that this curve becomes a limit cycle for a 
particular value of a. 

TuroremM lV. All cyeles of a given field containing the same critical points 
as a given Bendixson curve can be constructed by rotation operations of a suit- 
able complete family. 

1 wish to thank Professors S. Lefschetz and M. Schiffer for their interest 
and advice concerning this work. 


ON THE ADDITION AND MULTIPLICATION THEOREMS FOR 
SPECIAL FUNCTIONS* 
By C. TRUESDELL 
Greapvuare [nsrrrure sor INDIANA UNIVERSITY 


Communicated by T. Y. Thomas, October 16, 1950 


Seme years ago I called attention to the unification in the theory of 
special functions achieved through systematic exploitation of properties of 
the solutions of the F-equation! 


OF a) +1) 1) 
a: ot 


One virtue of the method is that many formal relationships are discovered 
and demonstrated in an automatic and trivial fashion, relationships which, 
if previously known at all, were originally derived at greater length and 
with more elaborate analytical tools. In particular, if Fiz + y, a) be an 
analytic function of y, we have by Taylor's theorem the addition formula 


Fiz, a + n), (2) 
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which contains as special cases most of the addition formulae and many of 
the generating expansions for familiar special functions.’ In this note I 


wish to point out that this same expansion may be written as a multiplica- 
tion theorem. First we put y — s for y in (2): 


>on (y 8)” 
Fly, a) =! is, + 29), (3) 


a n! 


If y = &z then (3) becomes the multiplication theorem 
— 1)"s" 
Fikz, a) = Fis, a + n), 


ae n. 


which is the desired result. 
If, in the general multiplication theorem (4), we put 


Fis, a) = (2s), (5) 
where /,(x) is the Bessel function of the first kind, after some reductions 
we readily obtain the known result* 

~~ (1 — /z\" 

J (ke) = k* (5) (6) 
n. « 
If J, be replaced by Y,,, (6) remains valid, If instead we put 
i 
F(s, a) = Ta + 1) (—2) LL (1/8) (7) 


where L,'°'(z) is the generalized Laguerre function, by substitution in (5) 
follows the first multiplication theorem of Erdélyi:* 


n=0\ 4 k 
while if we put 
Fiz, a) = (8), 
we similarly obtain the second multiplication theorem of Erdélyi:* 


ned 


(10) 


Of the at least eleven multiplication theorems which hold for the hyper- 
geometric functions, it may be of interest to write down two of the simplest, 
those resulting from the respective choices 


F(z, a) = F («. ), (11) 


\ 
| 
| 
s 
at (4) 
| 
- 
| 
\ 
[| i 
Ly 
| | 


74 MATHEMATICS: C. TRUESDELL Proc. N. A. 5. 


F(z, a) = (1 — s)'**~*F(b, c; a; 8), (12) 
namely, 
1 
Fla, bj ke) ‘) F(a + n, 6; c; s). 
(13) 
F(a, b; ¢ +m; 2). (14) 
Finally, we may put 
F(s, a) = I'(a)f(a, 2), (15) 


where {(a, 2) is the generalized Euler-Riemann zeta function, thus ob- 
taining the expansion 


t(a,ks) = > ( ) (1 — k)"s" + n, 2), (16) 
ae 

a result which when s = | yields a formula for {(a, &) as a series of simple 

zeta functions. 

We conclude with three remarks. The first is that since all the results 
of this note have been derived by simple substitutions in (2), without use 
of any limit operation, and since (2) is valid whenever the series on the 
right converges, it follows that all our results are valid subject only to 
convergence of the series involved. The second is that since & and z occur 
symmetrically on the left side of (4), by interchanging them upon the 
right side and equating the results we obtain the reciprocal theorem 


= Fiz,a +n) = Fik, a + n), (17) 
n=O n: 


subject to the convergence of both series. The third is that by putting 
k « Othe formula (4) enables us to find 


=~. (—)*s" 
ola) = FiO, a) = Fiz, a + n), (18) 


aw 
a result which is interesting both because the left-hand side is independent 
of sg and because the function ¢(a) occurs explicitly in many of the theorems 
concerning solutions of the F-equation. 


* Prepared under Navy Contract No. NSonr-180, Task Order V, with Indiana Uni- 
versity 
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* Truesdell, C., An Essay Toward a Unified Theory of Special Functions, Princeton, 
1948. 

2 Op. cit., 14. 

* Watson, J. N., A Treatise on the Theory of Bessel Functions, Cambridge, 1922. See 
45.22. 

Erdélyi, Funktionalrelationen mit konfluenten hypergeometrischen Funktionen,” 
1, Math. Z., 42, 125-143 (19386-1937), eq. (5,3). 

* Op. cit., eq. (6,2). In this note we employ the definition of the generalized La- 
guerre function given by E. Pinney, ‘‘Laguerre Functions in the Mathematical Founda- 
tions of the Electromagnetic Theory of the Paraboloidal Reflector,” J. Math. Phys., 
25, 49-79 (1946), so that the lower index is not restricted to integer values. Thus our 
results (8) and (10) are fully equivalent to Erdélyi’s formulae (5,1) and (6,1), which 
are expressed in terms of the M;, » functions of Whittaker. Generalizations of Erdélyi's 
formulae (5,4) and (6,3) may be obtained either by applying Kummer's first transforma- 
tion to our expansions (8) and (10) or by inserting successively the two solutions 


— a + 1) and — d) ) int the general multipli- 
3 


cation theorem (4). 
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